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EFFECT  OF  TRANSFAMINAR  REINFORCEMENTS  AND  HYBRIDIZATION  ON  DAMAGE  RESISTANCE 
AND  TOLERANCE  OF  COMPOSITE  LAMINATES 

ABSTRACT 

It  was  shown  that  the  damage  resistance  and  tolerance  of  laminated  composites 
can  be  enhanced  by  the  employment  of  translaminar  reinforcements  (TER)  such  as 
stitching,  z-pinning  and  3D  weaving  and  also  by  hybrid  composites. 

A  non-dimensional  analytical  model  focused  on  Mode  I  delamination  was 
developed  to  understand  the  role  of  the  TER  on  delamination  behavior.  An  explicit 
formula  for  the  apparent  interlaminar  fracture  toughness  was  derived  in  terms  of  the 
inherent  fracture  toughness  of  original  materials  and  the  bridging  force  due  to  z-pins. 

This  model  is  capable  of  estimating  the  apparent  fracture  toughness,  the  bridging  length 
and  allowable  bridging  force  thus  can  be  useful  in  the  design  of  TLR  for  composite 
laminates. 

Along  with  understanding  advantages  of  TLR  in  increasing  the  damage  tolerance 
of  laminated  composites,  the  damage  behavior  of  laminated  composites  subjected  to 
low  velocity  impact  loading  was  studied.  Based  on  the  similarity  in  damage  development 
between  quasi-static  and  dynamic  loadings  observed  through  the  short  beam  shear 
(SBS)  tests,  the  FE  analyses  of  the  SBS  specimens  for  quasi-static  indentation  and  at 
several  rates  of  low-velocity  impact  loadings  were  performed.  The  results  reveal  that 
inertia  effects  in  the  typical  velocity  range  of  the  striker  in  SHPB,  around  lOm/s,  can  be 
negligible,  and  hence  the  quasi-static  analysis  is  useful  and  valid  in  the  study  of  damage 
in  composite  specimens  under  low  velocity  impact  loading. 

The  delamination  behavior  of  3D  woven  composites  was  investigated  focusing  on 
the  effect  of  z-yam.  The  3D  woven  composites  containing  both  single  and  double  zl02 
yams  were  chosen  and  compared  with  the  2D  plain  woven  laminate.  The  double  z-yam 
woven  composite  exhibited  enhanced  damage  tolerance  compared  to  the  single  z-yam 
and  the  plain  woven  laminate.  The  relative  sliding  motion  between  two  layers  is 
constrained  by  z-yams  thus  the  crack  propagation  of  delamination  is  suppressed.  This 
mechanism  increases  the  apparent  interlaminar  fracture  toughness  of  the  composites. 

The  interlaminar  shear  stress  profiles  in  the  various  hybrid  composites  were 
obtained  from  linear  FE  analyses  and  compared  within  the  framework  of  damage 
resistance.  Some  hybrid  configurations  resulted  in  reduced  maximum  shear  stress 
value  for  a  given  contact  force  thus  demonstrating  higher  damage  resistance.  Although 
matrix  cracks  may  change  the  location  of  delamination  and  decrease  damage 
resistance,  hybrid  composites  can  still  be  superior  in  stmctural  applications. 
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EFFECT  OF  TRANSLAMINAR  REINFORCEMENTS  AND  HYBRIDIZATION  ON 
DAMAGE  RESISTANCE  AND  TOLERANCE  OF  COMPOSITE  LAMINATES 

By 

Min  Cheoi  Song 
December  2012 

Chair:  Bhavani  V.  Sankar 
Cochair:  Ghatu  Subhash 
Major:  Mechanicai  Engineering 

The  effects  of  translaminar  reinforcements  (TLR)  and  hybridization  on  impact 
damage  resistance  and  damage  toierance  of  iaminated  composites  were  analyzed.  The 
TLR  increase  the  apparent  fracture  toughness  of  composite  laminates  and  contribute  to 
improved  impact  properties.  Analyticai  and  computationai  methods  were  used  to 
investigate  the  damage  of  iaminated  composites  with  special  focus  on  delamination.  A 
non-dimensional  analytical  model  for  mode  I  delamination  of  z-pinned  composites  was 
developed  and  verified  using  finite  element  (FE)  analysis.  The  anaiytical  and  FE  modeis 
were  compared  with  experimental  results  to  evaluate  the  adequacy  of  the  model  in 
describing  the  role  of  translaminar  reinforcements.  The  influence  of  the  TLR  on  the 
apparent  fracture  toughness  and  bridging  length  was  quantitatively  investigated  through 
parametric  studies.  The  maximum  ailowable  bridging  force  before  inherent  faiiure  of  the 
materiai  was  suggested  as  weii. 

A  3D  woven  composite  was  chosen  to  study  the  effect  of  z-yarns  on  impact 
damage.  A  detaiied  analysis  was  performed  to  understand  the  iimitations  of  quasi-static 
anaiysis  in  low-veiocity  impact  studies.  The  static  equiiibrium  and  shear  stress  evolution 
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in  a  beam  made  of  orthotropic  material  subjected  to  quasi-static  and  impact  loadings 
under  various  impact  speeds  were  compared.  The  maximum  impact  velocity  for  which 
static  analyses  are  adequate  was  determined. 

Short  beam  shear  (SBS)  test  specimens  of  plain  woven  laminated  composite  and 
3D  woven  composites  were  analyzed  for  the  purpose  of  evaluating  the  effect  of  z-yarns 
on  delamination.  The  FE  model,  that  considered  both  intralaminar  and  interlaminar 
fracture,  predicted  damage  patterns  (transverse  cracks  and  delaminations)  observed  in 
the  tests.  The  z-yarns  in  the  composites  increase  damage  tolerance  by  interrupting 
crack  propagation. 

Finally,  the  advantage  of  hybridization  of  laminate  composites  was  evaluated 
using  the  FE  models.  Some  hybrid  composites  reduce  the  maximum  interlaminar  shear 
stress  in  beam-type  specimens,  and  therefore  enhance  the  damage  resistance  of 
composite  laminates.  Suggestions  for  continuing  this  study  in  the  future  are  included  at 
the  end. 
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CHAPTER  1 
INTRODUCTION 

Motivation 

Composite  materials  are  a  combination  of  two  or  more  materials  with  properties 
that  cannot  be  achieved  by  any  one  of  the  constituent  materials  [Gibson  1994]. 
Composite  laminates  have  been  wideiy  used  in  aerospace  and  automobile  structures, 
sporting  goods  and  military  equipment  because  of  their  superior  in-piane  stiffness  and 
strength.  However,  fiber  composites  are  anisotropic,  showing  reiativeiy  low  stiffness  and 
strength  in  directions  other  than  the  fiber  direction.  In  the  case  of  laminated  composites, 
delamination,  i.e.,  separation  of  layers  or  plies,  is  one  of  the  most  significant  damage 
modes,  and  is  a  key  element  of  consideration  in  the  design  and  operation  of  composite 
structures.  Toughened  resins  can  solve  this  problem  to  some  extent.  Rubber  [Kumar 
and  Singh  2000,  Yan  et  ai.  2002]  or  thermopiastic  toughened  epoxy  [Wong  et  al.  2010] 
is  an  example  of  such  an  approach.  Recently  nano-paticles  have  been  added  to  the 
resin  phase  to  improve  the  fracture  toughness  [Kalarikkai  et  al.  2006,  Kim  et  ai.  2008, 
BoesI  et  ai.  2011].  However,  the  improvement  in  interiaminar  fracture  toughness  due  to 
the  above  mechanisms  is  not  significant.  Another  method  of  improving  interiaminar 
properties  is  the  introduction  of  reinforcements  in  the  thickness  or  z-direction.  These  are 
called  translaminar  reinforcements  (TLR)  which  seem  to  increase  the  fracture 
toughness  and  strength  by  an  order  of  magnitude.  Thus,  there  seems  to  be  a 
tremendous  interest  in  TLR  for  laminated  composites. 

Stitching,  z-pinning  and  3D  weaving  are  some  of  the  methods  by  which  TLR  can 
be  provided  into  composite  iaminates.  Stitched  composites  can  be  manufactured  by 
passing  a  sewing  needle  with  a  thread  or  yarn,  looping  or  interlocking  with  bobbin 


13 


thread,  through  laminates  (Figure  1-1  A).  Several  researchers  have  studied  delamination 
resistance  and  damage  tolerance  of  stitched  composites  under  Mode  I,  Mode  II  and 
Mixed  Mode  loading  [Sharma  and  Sankar  1997,  Chen  et  al.  2005,  Dransfield  et  al. 

1998,  Sankar  and  Zhu  2000,  Jang  and  Sankar  2005,  Rys  et  al.  2010].  In  general, 
stitching  reduces  the  in-plane  stiffness  and  strength  properties  [Dransfield  1994,  Mouritz 
1997,  2004].  Although  stitching  was  found  to  reduce  the  impact  damage  resistance,  it 
significantly  improved  the  impact  damage  tolerance  [Sharma  and  Sankar  1997]  and 
compression  after  impact  (CAI)  [Farley  et  al.  1992]. 

Z-pinning  is  a  form  of  TLR  similar  to  stitching  in  which  rigid  rods  or  pins  are 
inserted  through  the  thickness  of  the  composite  structure  (Figure  1-1B).  It  is  an 
alternative  method  to  stitching  and  suitable  for  thick  composites  and  sandwich 
structures  [Wallace  et  al.  2001].  The  z-pinning  process  using  an  ultrasonic  insertion 
device  is  efficient  when  pre-pregs  are  used  [Lander  2008].  Z-pin  increased  the  ultimate 
strength  of  joints  [Byrd  and  Birman,  Chang  et  al.  2006]  while  degrading  in-plane 
stiffness  and  strength  [Mouritz  2007,  Steeves  2006,  Chang  et  al.  2007].  Increase  in 
interlaminar  fracture  toughness  in  z-pinned  composites  is  due  to  the  pull  out  mechanism 
of  the  pins  and  it  depends  only  on  the  frictional  and  cohesive  forces  between  the 
reinforcement  and  surrounding  matrix  material  [Cartie  2000,  Cartie  et  al.  2004,  2006, 

Dai  et  al.  2004].  However,  it  is  not  as  effective  as  stitching  as  the  ends  of  the  z-pins  are 
not  anchored  as  in  stitching. 

Although  stitching  is  a  practical  method  of  improving  interlaminar  properties,  the 
reduction  in  in-plane  properties  due  to  damage  to  the  composite  preforms  during  the 
stitching  process  has  become  a  matter  of  concern.  At  the  same  time  developments  in 
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textile  technology  has  lead  to  advanced  weaving  and  braiding  processes  which  can 
produce  three  dimensionally  reinforced  composites.  Textile  processes  are  well 
developed,  fast  and  easy  to  implement.  Complicated  performs  could  be  made  using 
textile  processes  where  separate  steps  to  implement  TLR  into  3D  woven  composites 
((Figure  1-1C)  are  not  necessary.  One  single  weaving  process  can  texture  warp  yarns 
(0-degree  tows),  weft  yarns  (90-degree  tows)  and  z  yarns  concurrently. 

Some  of  the  3D  performs  are:  3D  orthogonal  weave,  layer  to  layer  interlock,  and 
angle  interlock  (Figure  1-5).  From  early  on  it  has  been  recognized  that  3D  braided 
[Gong  and  Sankar  1991]  and  3D  woven  composites  [Brandt  et  al.  1996,  Byun  et  al. 
1990,  Walter  et  al.  2010]  improve  impact  properties  and  damage  tolerance  of 
composites.  However,  there  has  been  no  systematic  study  to  understand  the 
effectiveness  of  various  3D  architectures  mentioned  above  on  the  impact  properties  of 
composites.  This  can  be  studied  either  by  experimental  testing  -  by  testing  various 
woven  architectures  under  various  impact  loading  conditions  or  by  simulating  their 
behavior.  The  simulations  could  be  performed  using  either  analytical  models  or 
computational  models.  Analytical  models  are  based  on  several  simplifying  assumptions 
and  hence  tend  to  be  approximate.  However,  they  are  useful  in  revealing  the  important 
physical  mechanisms  that  are  responsible  for  the  improvements  in  the  impact  behavior. 
They  are  also  useful  in  optimization  studies  which  require  thousands  of  Monte  Carlo 
simulations.  Computational  models  use  the  finite  element  analysis.  They  tend  to  be 
more  realistic.  However,  they  are  expensive  and  time  consuming.  In  this  research  we 
propose  to  use  both  analytical  and  computational  models  to  understand  the  effect  of 
TLR  in  various  3D  woven  composites. 
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Hybrid  composites  also  are  capable  of  enhancing  impact  resistance  [Wang  et  al. 
1991]  and  interlaminar  fracture  toughness  [Hwang  and  Huang  1999].  Numerous 
combinations  with  a  variety  of  fibers  and  resins  are  possible  in  improving  composite 
properties  (Figure  1-6).  The  hybrid  composites  using  carbon  and  glass  woven  fabrics 
shows  superior  in-plane  properties  [Pandya  et  al.  2011]  and  CAI  [Naik  et  al.  2001] 
compared  to  glass  woven  fabric.  And  the  carbon/glass  hybrid  composites  have 
increased  failure  strain  [Stevanovic,  1992].  Combination  of  Kevlar  fiber  composites  and 
carbon  fiber  composites  possess  superior  ballistic  performance  in  armor  applications 
[Grujicic  et  al.  2006]. 

The  present  study  is  concerned  with  the  effect  of  TLR  on  the  damage  resistance 
and  damage  tolerance  of  composite  laminates.  Damage  resistance  (specifically  impact 
damage  resistance  in  the  present  context)  is  the  ability  of  the  material  to  function 
without  undergoing  significant  damage  as  exhibited  by  the  loss  of  stiffness  or  strength 
when  subjected  to  a  certain  load.  Damage  tolerance  pertains  to  material  that  is  already 
damaged.  Damage  tolerance  is  the  ability  of  the  material  that  has  sustained  some 
damage  to  continue  to  perform  and  carry  loads  albeit  with  reduced  stiffness  of  the 
structure.  The  advantages  of  hybrid  composites  -  both  inter-ply  and  intra-ply  hybrid 
composites  and  functionally  graded  composites  -  will  be  evaluated. 

Review  of  Research  Approach 

In  the  following  we  discuss  the  technical  challenges  In  solving  the  proposed 
problem  and  the  approach  we  will  take  to  address  those  challenges. 

The  major  Issues  in  the  analytical  model  of  TLR  composites  are  to  identify  the 
interactions  between  the  z-reinforcement  and  the  surrounding  composite  and  to 
represent  them  using  appropriate  mathematical  models.  In  computational  models, 
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damage  modes  should  be  defined  and  appropriate  numerical  strategy  to  represent  the 
damage  needs  to  be  developed.  The  way  to  model  delamination  between  layers  of 
composites  is  crucial  since  impact  resistance  and  damage  tolerance  are  significantly 
affected  by  delamination.  At  the  same  time  intralaminar  failure  should  be  considered 
since  this  failure  mode  seems  to  interact  with  interlaminar  delamination.  In  order  to 
characterize  the  effect  of  z-yarn  on  impact  resistance  and  damage  tolerance,  the 
mechanism  of  damage  initiation  and  evolution  will  be  investigated  by  simulation  of  the 
static  short  beam  shear  (SBS)  test,  since  static  indentation  tests  are  useful  in 
understanding  low  velocity  impact  behavior.  Experimental  studies  have  indicated  that 
damage  mechanisms  in  static  indentation  and  dynamic  impact  tests  are  similar  and  one 
can  learn  more  from  static  tests  about  impact  damage. 

Representation  of  Bridging  Zone 

One  of  the  biggest  characteristics  during  delamination  of  composites  reinforced  by 
TLR  is  the  development  of  the  bridging  zone  denoted  by  length  c  in  Figure  1-7. 

The  bridging  zone  can  be  regarded  as  a  process  zone  where  extra  energy  is 
dissipated  due  to  the  failure  of  TLR.  In  other  words,  the  energy  to  create  delamination  in 
transversely  reinforced  composites  is  higher  than  that  for  original  unreinforced 
composites.  This  increase  in  energy  required  increases  the  delamination  resistance. 

This  increase  due  to  TLR  can  be  quantified  by  comparing  fracture  toughness.  The 
relationship  between  the  two  values  of  the  fracture  toughness  -  that  of  reinforced 
composites  with  TLR  and  unreinforced  composites  -  is  given  by 

Gc-app=Gc+^  (1-1) 
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where  and  ,  respectively,  are  the  apparent  fracture  toughness  of  reinforced 

composites  with  TLR  and  the  fracture  toughness  of  the  original  composite;  and  ris  the 
energy  dissipation  rate  due  to  the  deformation  and  failure  of  TLR. 

So,  an  analysis  of  the  bridging  zone  can  provide  an  understanding  of  the  effect  of 
TLR  on  delamination.  The  behavior  of  the  entire  bridging  zone  must  be  characterized  by 
failure  mechanisms  or  a  bridging  law  for  each  TLR.  Experiments  [Tan  et  al.  2010,  Dai  et 
al.  2004]  have  been  conducted  to  measure  the  bridging  law.  It  has  been  found  that  a 
simple  force-displacement  relationship  of  the  z-reinforcement  could  describe  the  effect 
of  TLR  with  sufficient  accuracy.  Moreover,  bridging  laws  of  TLR  can  be  idealized  as 
linear  elastic  and  linear  softening  (linearly  decreasing  function)  laws  for  stitching  and  z- 
pinning,  respectively. 

The  use  of  spring  elements  is  very  adequate  and  efficient  since  their  behavior  can 
be  defined  by  a  force-displacement  relationship  and  makes  it  possible  to  simplify 
analysis.  The  spring  elements  with  these  bridging  laws  are  capable  of  being  directly 
implemented  into  double  cantilever  beam  specimens  (DCB)  or  end  notch  flexure  (ENF) 
specimens  since  the  relative  displacement  at  the  two  end  nodes  can  be  related  to  the 
deflection  of  beam.  So  the  DCB  and  ENF  specimens  incorporating  spring  elements 
have  been  used  for  representing  Mode  I  and  Mode  II  fracture  using  analytical  model  or 
finite  element  analysis  (FEA).  Discrete  spring  models  and  distributed  traction  models 
are  available  for  these  applications.  The  discrete  spring  models  [Ratcliffe  and  O’brien 
2004]  place  spring  elements  discretely  to  represent  TLR  while  the  distributed  traction 
model  [Sankar  and  Dharmapuri  1998,  Sankar  and  Zhu  2000,  Mabson  and  Deobald 
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2000,  Sridhar  et  al.  2002,  Robinson  and  Das  2004,  Byrd  and  Birman  2005]  uses 
continuous  elements  obtained  by  smearing  the  discrete  forces. 

Delamination  Modeling  using  FEM 

In  order  to  demonstrate  delamination  using  FEM,  J-integral  [Rice  1968]  and  virtual 
crack  closure  technique  (VCCT)  based  on  fracture  mechanics  [Rybicki  et  al.  1977, 
Krueger  2002,  GrassI  and  Zhang  2003]  or  the  cohesive  zone  method  (CZM)  based  on 
damage  mechanics  have  been  widely  used.  The  CZM  has  advantages  in  that  it  is  not 
necessary  to  define  a  pre-existing  crack  while  the  methods  based  on  fracture 
mechanics  require  pre-existing  cracks  with  which  strain  energy  release  rate  (SERR)  can 
be  calculated  and  used  in  a  fracture  criterion.  However,  The  CZM  is  also  related  to 
fracture  mechanics  since  dissipated  energy  of  cohesive  elements  is  the  same  as 
fracture  toughness  of  the  material,  and  this  fact  makes  it  possible  to  model 
delamination.  The  CZM  was  introduced  by  Barenblatt  [Barenblatt  1962]  and  Dugdale 
[Dugdale  I960].  The  traction-separation  law  is  a  commonly  adopted  approach  to 
describe  behaviors  of  the  cohesive  zone.  Several  traction-separation  laws  [Tvergaard 
and  Hutchinson  1992,  Xu  and  Needleman  1994,  Camacho  and  Ortiz  1996]  were 
suggested.  Among  them,  the  bi-linear  traction-separation  laws  [Geubelle  and  Baylor 
1998]  and  have  extensively  used  for  delamlnation  problems  in  Mode  I,  Mode  II  and 
Mixed  mode  problems.  The  bl-linear  traction-separation  law  Is  depicted  in  Figure  1-10 
and  can  be  defined  as 
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where  a  is  traction,  k  is  stiffness,  D  is  damage  variabie,  d  is  displacement,  is  the 
displacement  at  damage  initiation  and  d^  is  the  final  displacement  [Davila  et  al.  2007]. 

The  damage  behavior  of  the  cohesive  eiement  incorporating  the  bi-iinear  traction- 
separation  iaw  can  be  defined  by  parameters  -  cohesive  strength,  critical  displacement 
and  fracture  toughness  among  which  only  two  parameters  are  independent,  as  the  area 
under  the  stress-displacement  diagram  should  be  equal  to  the  fracture  toughness.  In 
addition  to  the  above  three  parameters,  the  initial  stiffness  or  the  penaity  stiffness  is 
also  a  crucial  parameter  that  affects  the  compliance  of  the  structure.  The  criteria  for 
damage  onset  and  propagation  should  be  defined  for  progressive  deiamination 
modeling.  Quadratic  stress  based  faiiure  criterion  [Cui  et  ai.  1992]  given  in  Eq.  (1-3) 
below  has  been  used  for  delamination  onset: 


(1-3) 


where  al  and  al  are  the  interlaminar  normai  and  shear  strengths. 

Power  law  criterion  [Wu  and  Reuter  1965]  (Eq.  1-4)  and  B-K  criterion  [Benzeggagh 
and  Kenane1996]  (Eq.  1-5)  are  weli  estabiished  for  damage  propagation 
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where  a  is  a  material  constant  and 


Gjc  -I-  {G,,^  Gjr ) 


A'' 


\Gj.  j 


=  G. 


(1-5) 


where  /;  is  a  material  constant  and  G^  =G,  +Gj, . 

Application  of  the  Bridging  and  Delamination  Model  for  a  Z-pinned  Composite 

The  FE  models  incorporating  the  cohesive  element  as  well  as  the  failure  of  TLR 
have  been  used  for  the  simulation  of  progressive  delamination  behavior  of  DCB 
reinforced  by  TLR  [Ratcliffe  and  Krueger  2006,  Dantuluri  et  al.  2007],  The  z-pinned 
composite  subjected  to  Mode  I  loading  has  been  studied,  since  the  debonding  of  skin- 
stiffener  T-joint  is  one  of  the  main  issues  for  designers.  This  anaiysis  can  provide  direct 
insights  in  understanding  the  role  of  the  TLR  to  damage  tolerance  by  evaluating 
apparent  fracture  toughness. 

A  summary  of  past  work  in  this  area  is  provided  in  Tabie  1  -1 .  In  the  following  we 
briefiy  describe  the  saiient  features  of  the  work  listed  in  that  table.  Cartie  [Cartie  2000] 
suggested  a  bi-linear  bridging  law  through  experiments  and  conducted  FE  simulation 
using  piain  strain  elements.  Robinson  and  Das  [Robinson  and  Das  2004],  Mabson  and 
Deobaid  [Mabson  and  Deobaid  2000]  and  Byrd  and  Birman  [Byrd  and  Birman  2005] 
used  a  linear  softening  bridging  model  whereas  Ratcliffe  and  O’Brien  [Ratciiffe  and 
O’Brien  2004]  proposed  a  discrete  element  analytical  model.  For  delamination  modeling 
of  z-pinned  DCB,  Dantuluri  et  al.  [Dantuluri  et  al.  2007]  used  2D  cohesive  eiements 
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incorporating  beam  elements  and  a  bi-linear  bridging  law.  Ratcliffe  and  Krueger 
[Ratcliffe  and  Krueger  2006]  used  3D  cohesive  model  with  solid  elements.  Grass!  and 
Zhang  [Grass!  and  Zhang  2003]  used  VCCT  for  delamination  modeling.  It  is  clear  from 
the  table  that  various  researchers  have  used  different  approaches  in  modeling  the 
stated  problem. 

Static  Indentation  and  Low-Velocity  Impact 

Numerous  studies  have  shown  [Kwon  and  Sankar  1993,  Nettles  and  Douglas 
2000]  that  static  indentation  tests  and  simulations  are  useful  in  understanding  low 
velocity  impact  behavior.  As  long  as  the  mass  of  the  impactor  is  higher  than  the  target 
mass,  beam  mass  in  the  present  case,  and  the  impact  velocity  is  much  less  than  the 
velocity  of  wave  propagation  in  the  target  medium,  the  impact  can  be  considered  as  low 
velocity  impact.  In  such  cases  static  tests  and  analyses  provide  good  information  on 
forces  at  which  damage  initiates  and  also  the  residual  stiffness  which  is  a  measure  of 
damage  tolerance. 

Damage  Behavior  in  Woven  Composites  during  Short  Beam  Shear  Test 

SBS  test  is  a  test  method  used  to  determine  interlaminar  shear  strength  (ILSS). 
SBS  tests  for  various  3D  woven  composites  were  conducted  to  investigate  the  effect  of 
z-yarns  on  ILSS  and  post-damage  behavior.  The  test  results  revealed  that  delamination 
is  a  predominant  damage  mode,  and  the  damage  tolerance  became  high  as  the  volume 
fraction  of  z-yarns  increased.  However,  with  the  increased  amount  of  z-yarns 
interlaminar  shear  strength  decreased  and  the  delaminations  were  contained,  leading  to 
post-damage  behavior  different  from  that  of  plain  woven  laminated  composites  [Walter 
et  al.  2010]. 
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One  of  the  significant  attributes  of  damage  evolution  was  that  matrix  cracks 
occurred  at  the  early  stage  of  loading  and  induced  interlaminar  delamination.  However, 
most  matrix  cracks  developed  inside  of  weft  yarns  (90-degree  tows)  as  shown  in  Figure 
1-12  where  ellipses  indicate  the  matrix  phase.  This  damage  pattern  is  a  common 
feature  in  cross-ply  laminated  composites  [0/90],  where  transverse  cracks  initiate  in  90- 
degree  plies  followed  by  delamination  between  90-degree  ply  and  0-degree  ply 
[Geubelle  and  Baylor  1998,  Ling  et  al.  2011].  In  fact  each  of  the  yarns  in  a  3D  woven 
composite  can  be  regarded  as  a  unidirectional  composite,  which  is  weak  in  the 
transverse  direction.  Therefore  the  effect  of  parameters  such  as  interlaminar  shear 
strength,  interlaminar  fracture  toughness  and  transverse  strength  of  yarns  on  damage 
initiation  and  evolution  should  be  investigated. 


Figure  1-1.  Schematics  of  composites  reinforced  by  the  TLR.  A)  Stitched.  B)  z-pinned. 
C)  3D  woven.  [Mouritz  2008]. 
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Figure  1-2.  NASA  and  Boeing’s  advanced  stitching  machine  [Lander  2008]. 


Figure  1-3.  Inserting  of  z-pin  [Cartie  2000,  Mouritz  2007]. 


Figure  1-4.  3D  weaving  machine  [Brandit  2001]. 
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Figure  1-5.  Schmatics  of  3D  woven  composites.  A)  3D  Orthogonal.  B)  Layer  to  layer 
interlock.  C)  Angle  interlock.  [Mahmood  et  al.  2011]. 


A  B 

Figure  1-6.  Hybrid  composites.  A)  Intraply.  B)  Interply. 


Figure  1-7.  Schematic  of  bridging  zone  developed  during  delamination 


Figure  1-8.  Bridging  laws  of  TLR.  A)  Linearly  increasing  (stitching).  B)  Linearly 
decreasing  (z-pinning). 
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A 
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Figure  1-9.  Spring  model  for  bridging  zone.  A)  discrete  spring  modei.  B)  distributed 
traction  model. 


Figure  1-10.  Traction-  separation  law  for  the  cohesive  element. 
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Figure  1-12.  Damages  after  SBS  test.  A)  In  plain  woven  laminated  composites.  B)  In  3D 
woven  composites.  [Walter  et  al.  2010]. 
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Table  1  -1 .  Summary  of  results  obtained  by  various  researchers  and  present  study. 


Reference 

Z-pin  in  Analytic  model 

Z-pin  model 
in  FEM 

Crack  propagation 
criterion  in  FEM 

Sub-laminate 
model 
in  FEM 

Comparison  with 
experimental  data 

Cartie  2000 

- 

Bi-linear  traction 

J  integral 

Plane  strain 

Yes 

Mabson  and  Deobald 
2000 

Distributed  force 
(Non-dimensional 
form-1) 

- 

- 

- 

- 

Byrd  and  Birman 

2005 

Distributed  force 
(Dimensional  form-1) 

- 

- 

- 

- 

Robinson  and  Das 
2004 

Distributed  force 
(Dimensional  form-2) 

- 

- 

- 

Yes 

Grass!  and  Zhang 
2003 

- 

Nonlinear  spring 
(Bi-linear  function) 

VCCT 

Shell  element 

Yes 

Ratcliffe  and  O’Brien 
2004 

Discrete  force 

- 

- 

- 

Yes 

Dantuluri  et  al.  2007 

- 

Nonlinear  spring 
(Bi-linear  function) 

2D  cohesive 
element 

Beam  element 

Yes 

Ratcliffe  and  Krueger 
2006 

- 

Nonlinear  spring 
(Bi-linear  function) 

3D  cohesive 
element 

Solid  element 

- 

Present  study 

Distributed  force 
(Non-dimensional 
Form-2) 

Nonlinear  spring 
(Linear  softening 
unction  and  tri-linear 
bridging  law) 

3D  cohesive 
element 

Shell  element 

Yes 
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CHAPTER  2 

ANALYSIS  OF  MODE  I  DELAMINATION  OF  Z-PINNED  COMPOSITE  USING  A  NON- 

DIMENSIONAL  ANALYTICAL  MODEL 

Our  goal  here  is  to  use  both  analytical  and  numerical  approaches  for  studying  the 
effect  of  z-pins  and  compare  them  with  avaiiabie  experimentai  data.  In  particular,  we 
have  developed  a  non-dimensional  model  that  will  be  useful  in  the  design  of 
transiaminar  reinforcements  for  a  given  appiication.  In  particular  we  have  used  a  3D 
shell  model  in  conjunction  with  3D  cohesive  eiements.  Both  linear  softening  laws  and 
more  realistic  tri-iinear  bridging  iaws  are  used  to  modei  the  z-pins.  In  order  to  clarify  the 
roie  of  z-pins,  a  non-dimensionai  anaiytical  modei  is  proposed.  The  solution  of  the  non- 
dimensionai  equation  uses  an  iterative  procedure.  We  have  also  derived  an  expression 
for  the  maximum  density  of  z-pins  that  can  be  ailowed  before  the  beam  faiis  otherwise. 
The  efficacy  of  the  anaiytical  model  is  verified  by  finite  eiement  simuiation  of  the  DCB 
specimen.  We  have  used  the  exampie  given  in  [Cartie  2000]  for  this  purpose.  In  the  FE 
simulation,  the  iigaments  of  the  DCB  are  modeled  using  shell  elements.  Cohesive 
elements  are  used  to  simulate  the  delamination  and  discrete  noniinear  elements  are 
used  to  model  the  z-pins.  The  agreement  between  the  analyticai  modei  and  FE 
simulations  is  found  to  be  exceiient  for  various  results  such  as  load-deflection,  load- 
crack  length  and  effective  fracture  toughness. 

It  was  found  that  steady-state  bridging  iength  and  maximum  apparent  fracture 
toughness  can  be  related  to  interlaminar  fracture  toughness  of  the  composite  and 
maximum  frictionai  force  supported  by  z-pins.  The  relationships  among  various 
parameters  are  non-dimensionalized  and  the  maximum  pin  friction  that  can  be  aiiowed 
before  the  composite  beam  itseif  fails  is  caicuiated.  The  non-dimensional  analytical 
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model  could  be  a  useful  design  tool  in  selecting  z-pins  for  composite  structures  to 
improve  interlaminar  fracture  toughness. 

Non-Dimensional  Anaiytical  Modei 

Consider  a  z-pinned  composite  DCB  specimen  of  thickness  2/7  with  initial  crack 
length  of  ao and  subjected  to  Mode  I  loading  as  shown  in  Figure  2-1  A.  A  pair  of 
transverse  forces  Fis  applied  at  the  tip  of  the  DCB.  Due  to  the  applied  transverse  force 
Fthe  initial  crack  tends  to  reach  the  current  crack  length  denoted  by  a  (Figure  2-1 B). 
When  the  crack  reaches  the  region  reinforced  by  z-pins,  a  bridging  zone  of  length  c 
begins  to  develop  (Figure  2-1 C).  We  can  define  an  apparent  crack  length  ap  which  is 
the  length  of  the  crack  up  to  the  beginning  of  the  bridging  zone.  Current  crack  length  (a) 
is  the  sum  of  the  apparent  crack  length  (ap)  and  the  bridging  length  (c).  When  the  z-pins 
start  to  be  completely  pulled  out  of  the  composite  as  shown  in  Figure  2-1 D,  the  bridging 
zone  becomes  fully  developed  and  a  new  fracture  surface  is  created  in  the  z-pinned 
zone.  Therefore  in  the  bridging  zone  the  pull-out  of  pins  is  partial.  Beyond  the  current 
crack  tip  the  pins  are  assumed  to  be  intact.  The  relative  deflection  at  the  tip  of  the  DCB 
is  denoted  by  Sas  shown  in  Figure  2-1 E.  Our  goal  is  to  determine  the  relationships 
among  F,  a,  cand  S.  Once  the  bridging  zone  length  cis  determined  other  parameters 
such  as  F,  a  and  S  can  be  found  by  solving  the  governing  equations.  In  the  following  we 
describe  the  non-dimensional  equation  governing  the  bridging  zone  and  procedures  to 
determine  the  bridging  length  for  a  given  load. 

We  assume  that  the  shear  deformation  is  negligible  and  use  Euler-Bernoulli  beam 
equations  to  model  the  ligaments  of  the  DCB  as  shown  in  Figure  2-2B.  We  also  assume 
that  z-pins  are  rigid  and  there  is  sufficient  friction  between  the  z-pins  and  the  composite 
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material  surrounding  it.  For  simplicity,  the  relation  between  the  frictional  force  {f)  and 
pull  out  or  slip  distance  (c/s)  is  idealized  as  a  linear  softening  function  as  shown  in  Figure 
2-3.  The  validity  of  this  assumption  will  be  later  verified  in  FE  simulation.  When  the  pins 
are  intact  they  can  exert  a  maximum  friction  force  of  fm-  As  the  pins  pull  out  of  the 
material,  the  loss  of  friction  Is  proportional  to  the  pullout  distance  ds-  When  the  pin  is 
completely  pulled  out  of  the  beam,  the  friction  force  reduces  to  zero.  Thus  the  f-ds 
relationship  is  given  by 

V  h  J 

where  h  is  half  the  thickness  of  DCB. 

Although  the  resistance  offered  by  the  pins  on  the  beam  is  discreet,  for  the 
purpose  of  the  analytical  model  we  smear  the  discrete  pin  resistance  as  continuous 
distributed  traction  (p)  acting  on  the  crack  surfaces  as  shown  in  Figure  2-2C.  Then  the 
traction  can  be  derived  as 

f  d\  (  d\ 

p  =  Nf  =  Nf^  l—f-  =p^  1-^  (2-2) 

y  hj  y  hj 

where  /Vis  z-pin  density  expressed  as  number  of  z-pins  per  unit  area. 

The  Euler-Bernoulli  beam  equation  for  one  of  the  ligaments,  say  upper  ligament 
(Figure  2-2C),  of  the  DCB  can  be  written  as: 

EI^  =  -bp  {2-3) 

dx 

where  b  is  the  beam  width.  The  effective  bending  rigidity  of  one  of  the  ligaments  of  the 
DCB  is  represented  by  the  term  El.  Note  that  the  flexural  rigidity  for  a  laminated 
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composite  can  be  taken  from  the  bending  stiffness  matrix  in  classical  plate  laminate 
theory.  Substituting  for  pfrom  Eq.  (2-2)  into  Eq.  (2-3),  we  obtain 


El 


-bp^^  =  -bp„ 
h 


(2-4) 


The  pullout  length  Is  equal  to  2w,  where  w  Is  the  deflection  of  the  top  or  bottom 


beam.  Hence,  the  governing  equation  takes  the  form 


El 


W 

■^PmT  =  -bPn 
h 


(2-5) 


The  origin  of  the  x-coordinate  is  assumed  to  be  the  point  at  which  the  bridging 
zone  begins.  One  should  note  that  the  origin  moves  as  the  crack  propagates.  The  four 
boundary  conditions  (BCs)  for  the  bridging  zone  are: 


w(0)  =  Wq,(0 <  Wq  < /z/ 2) 

d^w 


w(c)  =  0 

-W=o 

dx^  ’ 


dx 


(0)  =  a,^(0) 

dx 


(2-6) 


In  the  above  equation  wq  is  a  prescribed  deflection  at  the  origin  of  the  coordinate 
system  as  shown  in  Figure  2-2C.  This  value  will  be  less  than  h/2  in  the  beginning  and 
will  increase  to  a  maximum  value  of  h/2  as  the  DCB  is  loaded.  Note  that  the  slip 
distance  of  the  z-pin  is  equal  to  the  total  opening  of  the  DCB  at  that  location,  i.e. 
2w(x)=ds(x).  As  the  crack  propagates,  the  bridging  zone  will  also  move  with  the  crack, 
but  the  crack  opening  displacement  will  remain  as  h  at  the  origin  with  2w(0)=ds(0)=h. 
The  terms  Van6  M,  respectively,  are  the  bending  moment  and  transverse  shear  force 
on  the  beam  cross  section. 
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One  should  note  that  the  bridging  length  c  is  still  an  unknown.  It  can  be 
determined  from  the  fact  that  the  strain  energy  release  rate  at  the  actual  crack  tip 
should  be  equal  to  the  Mode  I  fracture  toughness  at  the  instant  of  crack  propagation. 
The  energy  release  rate  can  be  determined  from  the  equation  derived  by  Sankar  and 
Sonik  [Sankar  and  Sonik  1995]  for  beam-like  specimens  containing  delamination. 


_  El  (dMc)'" 
bEI  b  dx^  ^ 


(2-7) 


Thus  the  condition  for  determining  c  is 


d^w(c)^ 
b  dx^  j 


(2-8) 


Before  we  solve  the  above  equations  we  will  non-dimensionalize  the  equations 
and  BCs  appropriately.  Normalizing  the  length  dimensions  by  h  and  forces  by  Elf,  the 
governing  equation  and  the  BCs  take  the  following  form: 


d'^w 

djf 


^Pj^  =  -P. 


(2-9) 


wherejc  =  -,a  =^,c  =  -,w  =  -  and 

h  ^  h  h  h  E 

The  equation  for  determining  c  (Eq.  (2-8))  takes  the  form 


^  d^w{c)^ 

^  dx^  j 


(2-11) 
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^  12G 

where  the  non-dimensional  fracture  toughness  is  given  by  = - — 

Eh 

The  solution  for  the  governing  Eq.  (2-9)  is 

vi>(x)  =  Cj  cos  Xx  +  C2  sin  Ax  +  C3  cosh  Ax  +  sinh  +  ~  (2-12) 

where  ^  ~ 

The  boundary  condition  at  the  point  x  =  0  varies  since  deflection  at  this  point 

increases  from  zero  at  the  beginning  of  loading  to  0.5  when  the  bridging  zone  is 
completely  developed.  Once  the  bridging  zone  is  completely  developed  this  value 

remains  constant  at  0.5  with  additional  increment  of  This  is  because  the  bridging 
zone  is  fully  developed  and  It  moves  with  the  crack  tip  as  it  advances. 

The  procedures  to  solve  the  above  set  of  equations  are  shown  in  the  flow  chart 
depicted  in  Figure  2-4.  The  initial  data  Includes  the  beam  properties,  characteristics  of 
the  z-pins  and  the  fracture  toughness  G/c.  The  deflection  at  the  beginning  of  the 
bridging  zone  (x  =  0)  begins  to  increase  as  the  load  is  applied. 

When  vv(0)=  0.5,  the  bridging  zone  is  fully  developed  and  hereafter  the  deflection 
of  the  beginning  of  the  bridging  zone  is  constant.  However,  the  apparent  crack  length 
increases  during  crack  propagation.  We  need  to  use  an  iterative  procedure  as  the 
bridging  length  (c)  is  not  known  a  priori.  The  strain  energy  release  rate  condition  at  the 
right  end  of  the  bridging  zone  ( x  =  c )  as  given  by  Eq.  (2-1 1 )  is  then  used  to  check  for 
correct  value  of  c  .When  a  given  c  satisfies  Eq.  (2-11),  then  the  procedure  to  determine 
c  is  terminated.  After  bridging  zone  is  fully  developed,  bridging  length  corresponding  to 
every  increment  of  crack  length  can  be  determined. 
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Verification 


A  Sample  Problem  for  Verifying  the  Non-Dimensional  Analytical  Model 

In  order  to  verify  the  analytical  model  a  z-pinned  composites  designed  by  Cartie 
[Cartie  2000]  as  shown  in  Figure  2-5  was  selected.  This  z-pinned  composite  was  also 
used  in  many  previous  works  [Cartie  2000,  Grassi  et  al.  2003].  The  configurations  and 
material  properties  are  represented  in  Table  2-1 .  These  values  were  used  for  both  the 
analytical  model  and  FE  simulation  presented  in  the  next  section. 

First,  the  load-deflection  curve  (Figure  2-6)  was  found  for  ap  values  ranging  from 
52.25  mm  to  62.25  mm.  Meanwhile  the  crack  length  as  the  summation  of  the  apparent 
crack  length  (ap)  and  the  bridging  length  (c)  obtained  from  the  procedure  depicted  in 
Figure  2-4  was  computed.  The  variation  of  the  bridging  length  as  a  function  of  DCB 
deflection  ^is  shown  in  Figure  2-7.  The  bridging  length  initially  increases  with  loading 
until  the  bridging  zone  is  fully  developed,  where  the  bridging  length  has  the  maximum 
value.  In  the  beginning,  the  apparent  crack  length  remains  constant  and  the  crack 
propagation  is  only  due  to  evolution  of  the  bridging  zone.  Once  the  bridging  zone  is  fully 
developed  there  is  little  change  in  the  bridging  length  thus  the  crack  propagation  is 
almost  due  to  the  apparent  crack  length.  In  fact  there  is  a  slight  decrease  in  the  bridging 
length  after  it  reaches  the  maximum  as  the  bending  moment  due  to  the  transverse  loads 
applied  is  a  function  of  the  apparent  crack  length.  If  a  pair  of  couples  is  applied  instead 
of  a  pair  of  forces  as  in  the  standard  DCB  test,  then  the  bending  moment  will  remain 
constant  as  the  crack  propagates,  and  one  can  see  a  steady  state  bridging  length. 

The  force  Facting  on  the  DCB  and  opening  displacement  SaX  the  end  of  the  beam 
can  be  obtained  using  the  following  relations: 
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II 

1 

II 

O 

(2-13) 

„  J  dw  1  d^w  ^ 

5  =  2  w{Qi)  +  a  +  a\ 

^  dx  ^  dx  ^ 

V  ^=0  x=Q  J 

(2-14) 

The  variations  of  transverse  force  and  crack  length  with  increments  in  opening 
displacement  are  shown  in  Figure  2-6.  Initially  the  transverse  force  increases  with 
opening  displacement  and  begins  to  drop  as  interlaminar  crack  in  unreinforced  region 
propagates.  However,  the  transverse  force  increases  during  the  development  of  the 
bridging  zone  and  decreases  again  with  the  movement  of  the  fully  developed  bridging 
zone  and  the  new  crack  surface. 

In  order  to  estimate  the  apparent  fracture  toughness  during  crack  propagation  we 
use  the  bending  moment  at  the  apparent  crack  tip  as  a  measure  of  increased  fracture 
toughness.  This  is  similar  to  calculating  the  fracture  toughness  of  the  unreinforced  beam 
using  Eq.  (2-7).  Then  the  apparent  fracture  toughness  is  defined  as 


_El(d^w{Qi)^ 

IC-app--^y  ^ 


(2-15) 


The  above  relation  can  be  non-dimensionalized  as  : 


^"w(0)^ 
dx^  ^ 


(2-16) 


The  variation  of  apparent  fracture  toughness  during  crack  propagation  is  shown  in 
Figure  2-8.  The  apparent  fracture  toughness  increases  during  development  of  the 
bridging  length  since  the  bending  moment  required  at  the  beginning  of  bridging  zone  to 
overcome  the  bridging  force  increases  gradually. 
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Maximum  attainable  apparent  fracture  toughness  can  be  predicted  by  the 
relationship  expressed  in  terms  of  interlaminar  fracture  toughness  and  bridging  force 
due  to  z-pins.  The  relationship  is  exactly  the  same  as  the  equation  based  on  energy 
balance.  As  the  crack  propagates  it  has  to  overcome  the  frictional  forces  in  the  z-pins. 
The  amount  of  extra  work  done  is  equal  to  the  area  under  the  load-deflection  diagram  in 
Figure  2-3.  Thus  the  maximum  apparent  fracture  toughness  can  be  derived  as: 

G,e-.„.,^=G„+N[\f^h\  (2-17) 

J 

Multiplying  throughout  by  (I2/E/7)  we  obtain  the  above  relation  in  a  non- 
dimensional  form  as: 

^IC-app-max  ~  ^Ic 

This  maximum  value  can  be  realized  when  a  pair  of  couples  is  applied  instead  of  a 
pair  of  transverse  forces.  As  shown  in  Figure  2-8,  the  apparent  fracture  toughness  for 
transverse  loading  is  always  less  than  the  maximum  value.  Furthermore,  our  definition 
of  apparent  fracture  toughness  in  Eq.  (2-16)  is  also  slightly  different  from  the  energy 
based  relation  given  in  Eq.  (2-18)  making  the  apparent  fracture  toughness  under 
transverse  loading  less  than  the  maximum  realizable  value.  This  is  useful  in  evaluating 
the  increase  in  apparent  fracture  toughness  for  a  given  design  of  z-pins.  Moreover  this 
value  is  related  to  strain  exerted  in  the  ligaments  of  DCB  and  will  be  discussed  later. 

A  FE  Simulation  for  Verifying  the  Non-Dimensional  Analytical  Model 

For  the  sake  of  comparison  with  the  analytical  model,  FE  simulation  of  the  same 
specimen  using  the  finite  element  software,  ABAQUS®,  was  performed.  The  FE  model 
containing  discrete  z-pins  can  verify  the  analytical  model  where  z-pins  are  smeared  and 
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represented  by  distributed  traction.  The  three-dimensional  FE  analysis  was  used  to 
simulate  crack  propagation  in  the  DCB  specimen  with  z-pins.  The  specimen  was 
modeled  using  shell  elements  (S4),  and  cohesive  elements  (COH3D8)  were  used  to 
simulate  progressive  delamination. 

The  behavior  of  cohesive  elements  can  be  characterized  by  a  bi-linear  traction 
separation  law  (Figure  2-9)  [Davila  et  al.  2007]  given  by 


a  =  {\-D)Kd 

where  (2-19) 


D  =  \ 


0 

i^d-do) 
d {dj^  -do^ 


,d  Kd^ 

,dQ<d  <dj- 


1  ,d^  <d 

where  cr is  traction,  K\s  stiffness,  D  is  damage  variable,  dls  displacement,  do  is 
displacement  at  damage  initiation  and  df  is  final  displacement. 

The  above  parameters  were  taken  as:  K=^Q^  N/mm,  crm=35Mpa  and  G/c=0.258 
N/mm  [Cartie  2000,  Davila  et  al.  2007]. 

The  nonlinear  spring  element  (CONN3D2),  whose  behavior  is  defined  by  a  linearly 
decreasing  force  (Figure  2-3)  for  each  z-pin,  was  also  implemented  between  the  two 
ligaments  of  the  DCB  (Figure  2-10). 

The  load-deflection  curves  and  the  variations  of  crack  lengths  with  deflection  are 
shown  in  Figure  2-6.  Both  the  FE  and  analytical  model  results  are  presented  with  that 
from  the  experiment  by  Cartie  [Cartie  2000].  The  agreement  between  the  analytical 
model  and  the  FE  simulations  is  satisfactory  for  both  load-deflection  and  delamination 
length.  The  slight  discrepancy  in  the  initial  slope  of  the  load-deflection  curve  between 
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the  analytical  model  and  the  FE  model  is  due  to  the  fact  that  the  analytical  model  uses 
Euler-Bernoulli  beam  theory  whereas  the  FE  model  considered  shear  deformation  as  it 
occurs  in  test  specimens.  The  good  comparison  between  the  analytical  model  and  the 
results  from  FEA  and  experiments  suggests  that  the  discrete  bridging  force  can  be 
represented  by  a  distributed  traction. 

Discussion  of  Results  from  the  Non-Dimensional  Model 
Parametric  Studies  using  the  Non-Dimensional  Analytical  Model 

Our  goal  was  to  study  the  effects  of  inherent  interlaminar  fracture  toughness  of  the 

^  h 

composite  material  ^‘^and  the  non-dimensional  frictional  force  on  (i)  the  maximum 

G 

apparent  fracture  toughness  ic-app-mwc steady  state  bridging  length  ^ .  Such 
relationships  are  extremely  useful  in  the  design  process  to  evaluate  the  influence  of 
design  variables  on  performance  of  the  composite.  In  this  parametric  study  we  assume 
that  the  DCB  is  loaded  by  end  couples  instead  of  transverse  forces.  This  assures 
steady  state  crack  propagation  in  the  beam  and  the  effect  of  increasing  crack  length  on 

the  results  is  thus  eliminated.  In  the  simulations,  *  is  varied  from  1 0'^  to1 and 
ranged  from  10'^°  (almost  zero  friction  representing  the  case  of  an  unreinforced  beam) 
to  0.1. 

Figure  2-1 1  shows  that  the  steady  state  bridging  length  decreases  with  increasing 

*  and  .  If  the  composite  material  is  inherently  tough,  it  will  not  allow  a  larger 
bridging  length.  Similarly  a  large  friction  force  also  will  reduce  the  bridging  length.  The 

Q 

maximum  apparent  fracture  toughness  (  js  another  measure  to  evaluate  the 

effect  of  translaminar  reinforcement.  Figure  2-12  shows  that  the  maximum  apparent 
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fracture  toughness,  computed  by  Eq.  (2-16),  varies  linearly  with  increasing 


and  P" 


This  result  is  in  accordance  with  Eq.  (2-18).  Since  the  range  of  is  wider  than  that  of 


G  G  ^ 

in  reality,  is  more  sensitive  to  the  value  of  .  In  other  words  the  apparent 

fracture  toughness  is  dominated  by  translaminar  reinforcement.  From  the  results  shown 
in  Figures  2-1 1  and  2-12,  it  is  clear  that  higher  friction  force  exerted  by  the  z-pins 
provides  enhanced  fracture  toughness  and  at  the  same  time  reduces  the  bridging 
length.  The  latter  is  important  to  maintain  the  stiffness  of  the  structure  for  a  larger 
bridging  length  leads  to  reduction  in  the  stiffness  of  the  structure. 

Maximum  Allowable  Translaminar  Reinforcement 

Even  though  large  frictional  force  between  the  z-pin  and  the  surrounding  matrix 
material  is  desirable  for  increased  fracture  toughness,  a  frictional  force  beyond  a  critical 
value  will  cause  the  beam  to  fail.  The  maximum  normal  strain  in  a  beam  cross  section  is 
given  by 


s 


max 


h  d^w 
2  dx^ 


(2-20) 


Note  the  strain  is  already  non-dimensional  and  the  right  hand  side  of  the  above 
equation  can  be  written  as 


s 


max 


1 

2 


d^w 

dx^ 


K 

~2 


(2-21) 


where  ^is  the  non-dimensional  curvature  which  has  the  maximum  value  at  x  =  0.  Using 
Eqs.  (2-16)  and  (2-18)  we  obtain 


k  = 


(2-22) 
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Let  us  assume  the  allowable  strain  in  the  composite  is  given  by 
Then 

^4  + ip,.  < 2p.  (2-23) 

From  the  above  equation  one  can  derive 

P.,<2K"-G„) 

or  (2-24) 

3  h 

The  above  equation  provides  an  upper  limit  on  the  z-pin  density  which  should  be 
taken  into  consideration  in  the  design  of  translaminar  reinforcements.  Thus  the 
allowable  density  or  maximum  frictional  force  of  z-pins  can  be  determined  at  a  given 
geometry  and  material  properties  of  a  composite  material. 

In  order  to  verify  the  result  for  the  maximum  allowable  translaminar  reinforcement 
given  in  Eq.  (2-24),  we  performed  finite  element  analysis  of  a  specimen.  The  properties 

used  were:  G,^  =0.258N/mm,/r=1.6mm,£'  =  138GPa,A(f^  =1.58  N/mm^[Cartie  2000]. 

The  FE  simulation  was  used  to  calculate  the  maximum  strain  in  the  ligaments  of 
the  DCB  specimen  at  the  instant  of  crack  propagation,  and  it  was  compared  with  the  eu 
obtained  from  Eq.(2-24).  Two  loading  cases,  namely  end  couples  and  transverse 
forces,  were  considered.  For  the  set  of  properties  used  Eq.(2-24)  yields  a  maximum 

strain  Su  =  4.5x10'^.  Note  that  if  Nf^  >1.58  N/mm^  then  the  maximum  strain  exceeds  the 

ultimate  strain.  From  the  FE  simulations  we  obtained  £max  =4.521  x  10'^  for  the  case  of 
end  couples  and  Smax  =3.906x10'^  for  the  case  of  transverse  loading  of  the  DCB  (Figure 
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2-13).  The  results  indicate  the  strains  in  the  beam  are  higher  for  the  moment  loading 
compared  to  that  due  to  the  transverse  loading.  However,  in  a  practical  structure  the 
delamination  will  experience  a  combination  of  shear  and  moment  loadings,  and  hence 
the  conservative  value  should  be  used.  That  means  the  maximum  allowable  stitch 
density  given  in  Eq.  (2-24)  should  be  used  although  it  is  applicable  only  to  moment 
loading  case,  and  hence  conservative. 

Effect  of  Tri-Linear  Bridging  Law 

So  far  we  have  used  the  linear  softening  bridging  law  for  developing  the  non- 
dimensional  analytical  model  and  verifying  the  same  using  FEA.  Although  this  bridging 
law  is  simple,  it  may  not  be  realistic.  Dai  et  al.  [Dai  et  al.  2004]  performed  pull-out  tests 
using  z-pins  of  various  diameters  to  determine  the  actual  bridging  law.  They  found  that 
a  high  value  of  debonding  force  was  reached  before  the  debonding  of  the  pins  began. 
After  the  debonding  was  initiated,  the  pull-out  force  dropped  to  a  lower  value  before 
reducing  linearly  to  zero  value  as  the  pins  were  pulled  out  steadily  against  frictional 
force  that  seemed  to  have  a  constant  coefficient  of  friction.  They  represented  this 
pullout  behavior  by  a  tri-linear  bridging  law  as  depicted  in  Figure  2-14.  The  tri-linear 
behavior  can  be  attributed  to  an  elastic  deformation  (from  (0,  0)  to  (0.0185,  35.3)  in 
Figure  2-14),  debonding  from  surrounding  matrix  (from  (0.0185,  35.3)  to  (0.17, 14.86)  in 
Figure  2-14)  and  slip-out  of  z-pins  (from  (0.17,  14.86)  to  (1 .6,  0)  in  Figure  2-14). 

The  finite  element  simulations  of  the  DCB  described  above  were  repeated  with  the 
tri-linear  bridging  law.  We  did  not  an  attempt  analytical  solution  as  it  will  be  very 
complicated  due  to  the  piecewise  continuous  bridging  law.  The  numerical  values  for  the 
tri-linear  law  were  chosen  such  that  the  area  under  the  force-displacement  curve,  which 
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is  related  to  the  increase  in  fracture  toughness,  is  equal  to  the  linear  softening  law  used 
in  the  earlier  example. 

The  areas  under  the  two  force-displacement  curves  in  Figure  2-14  are  equal  to 
each  other.  The  resulting  load-deflection  diagram  of  the  z-pinned  DCB  specimen  is 
shown  in  Figure  2-15.  The  load  deflection  behavior  is  almost  identical  to  that  obtained 
using  the  linear  bridging  law  indicating  that  the  apparent  fracture  toughness  is  same  in 
both  cases.  It  must  be  noted  that  the  peak  force  in  the  tri-linear  case  is  almost  twice  as 
that  of  the  linear  softening  law.  The  variations  of  crack  length  as  a  function  of  the  DCB 
specimen  opening  (deflection)  are  also  similar  in  both  cases.  This  indicates  that  the 
details  of  the  bridging  law  do  not  seem  to  affect  the  global  behavior  of  the  specimen  as 
long  as  the  energy  dissipated  by  the  pins  is  properly  accounted  for. 

Summary  and  Conclusions 

Mode  I  delamlnation  propagation  in  DCB  specimens  containing  z-pins  is  studied.  A 
simple  analytical  model  based  on  a  linear  softening  type  bridging  law  for  the  z-pIns  has 
been  developed  and  suitable  non-dimensional  parameters  have  been  Identified.  The 
load-deflection  curve  of  the  DCB  specimen  was  calculated  using  the  analytical  model.  It 
is  seen  that  the  bridging  zone,  wherein  the  pins  are  partially  pulled  out,  develops  as  the 
crack  propagates,  but  attains  a  steady  state  value.  The  length  of  the  bridging  zone  is  a 
function  of  the  Mode  I  fracture  toughness  and  the  frictional  force  between  the  z-pIns  and 
the  surrounding  material.  An  expression  was  derived  for  the  apparent  or  effective 
fracture  toughness  values.  Although  increase  in  frictional  force  as  the  z-pins  increases 
the  fracture  toughness,  there  is  an  upper  limit  to  this  friction  as  the  DCB  ligaments 
would  break  If  the  friction  Is  very  high.  The  limiting  value  of  the  pin  friction  is  derived. 
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The  efficacy  of  the  analytical  model  was  evaluated  by  the  simulation  of  the  DCB 
specimen  using  finite  element  simulations.  In  the  FE  model  the  delamination 
propagation  was  simulated  by  cohesive  elements  and  the  z-pins  were  modeled  as 
discrete  nonlinear  elements.  The  results  for  the  load-deflection  curve  and  the  crack 
bridging  zone  length  agreed  quite  well  with  the  analytical  model.  As  an  alternative  to  the 
linear  bridging  law,  a  more  realistic  tri-linear  bridging  law  was  used  in  the  FE 
simulations.  It  is  found  that  the  global  delamination  behavior  of  the  specimen  was  not 
affected  much  as  long  as  the  energy  dissipated  by  the  pins  is  kept  the  same. 

The  non-dimensional  model  with  few  parameters  will  serve  as  a  design  tool  when 
translaminar  reinforcements  such  as  z-pins  are  selected  for  laminated  composite 
structures  in  order  to  improve  their  fracture  toughness.  The  analytical  models  will  also 
be  useful  in  optimization  studies  and  simulation  of  large  composite  structures  containing 
translaminar  reinforcements. 
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Figure  2-1 .  Schematic  of  initial  geometry  and  development  of  bridging  zone  by  z-pins  in 
the  composite  DCB. 
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Figure  2-2.  Idealization  of  bridging  zone  using  a  beam  model. 
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Figure  2-3.  Force-displacement  relation  of  the  z-pin. 


Figure  2-4.  Flowchart  of  the  procedures  for  solving  the  non-dimensional  governing 
equation. 
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Nondimensional  load 


Figure  2-5.  DCB  specimen  reinforced  by  z-pins  [Cartie  2000]. 
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Figure  2-6.  Load  and  crack  length  variation  as  a  function  of  DCB  deflection. 


47 


Nondimensional  crack  length 


Figure  2.-1 .  Bridging  iength  and  crack  iength  as  a  function  of  DCB  deflection. 


Figure  2-8.  Variation  of  apparent  fracture  toughness  during  crack  propagation. 
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Figure  2-9.  Traction-  separation  iaw  for  the  cohesive  eiement. 
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Cohesive  element 

Figure  2-10.  Cohesive  and  spring  elements  in  the  FE  model  of  the  DCB. 


Figure  2-1 1 .  Non-dimensional  steady  state  bridging  length. 
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Figure  2-12.  Apparent  fracture  toughness  as  a  function  of  maximum  friction  force  and 
interlaminar  fracture  toughness. 
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Figure  2-13.  Comparison  of  strain  with  different  ioading  types. 


Figure  2-14.  The  tri-linear  bridging  iaw  is  indicated  by  the  solid  line.  The  dotted  line  is 

the  iinear  softening  law  in  Figure  2-3.  The  areas  under  the  force-dispiacement 
diagram  for  both  laws  are  the  same. 
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Figure  2-15.  Load-deflection  curve  from  various  bridging  laws. 

Table  2-1 .  The  various  dimensions  and  properties  of  the  DCB  used  in  the  numerical 


simulation  [Cartie  2000] 


Dimension  or  Property 

Value 

B 

20  mm 

H 

1.6  mm 

E, 

138  Gpa 

E2 

11  GPa 

Vl2 

0.34 

Gi2 

4.4  GPa 

Gic 

258  N/m 

do 

49  mm 

Fm 

18.43  N 

z-pin  density 

0.5  % 

z-pin  diameter 

0.28  mm 
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Crack  length  (mm) 


CHAPTER  3 

APPLICABILITY  OF  QUASI-STATIC  ANALYSIS  FOR  LOW  VELOCITY  IMPACT 
Short  Beam  Shear  (SBS)  tests  are  widely  used  to  characterize  and  quantify  impact 
resistance  of  the  composites  with  various  constituent  materials  and  fiber  architectures 
for  both  quasi-static  and  dynamic  loading.  THE  SBS  test  is  basically  a  three-point  bend 
test  with  a  much  shorter  span  length.  The  ratio  of  length  to  height  of  the  specimen  {L/h 
ratio)  is  typically  less  than  10.  The  test  designated  as  ASTM  D2344  was  originally 
designed  for  measuring  interlaminar  shear  strength  of  laminated  composites.  Since  it 
closely  mimics  low-velocity  impact  response  and  damage,  it  has  been  used  to  measure 
impact  resistance  and  impact  damage  tolerance  of  composite  materials.  The 
interlaminar  shear  strength  (ILSS)  is  calculated  from  the  maximum  force  Pmaxbefore 
failure  initiates  in  the  specimen. 

The  interlaminar  shear  strength  is  calculated  as 

ILSS  =  -^  (3-1) 

4  bh 

The  above  equation  is  based  on  parabolic  variation  of  transverse  shear  stress  through 
the  thickness  of  the  beam.  It  is  true  only  for  homogeneous  materials  and  also  the  cross 
section  should  be  away  from  the  loading  region.  The  above  equation  which  is  based  on 
static  loading  is  also  being  used  for  understanding  impact  damage.  It  assumes  that  the 
specimen  is  in  quasi-static  equilibrium,  and  neglects  any  dynamic  effects.  However,  the 
above  assumption  would  be  valid  only  below  certain  impact  velocity.  Hence,  we  perform 
finite  element  simulations  of  dynamic  loading  cases  and  calculate  the  detailed  stress 
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field  and  compare  the  results  with  corresponding  static  solution  in  order  to  verify  the 
validity  of  quasi-static  assumption  for  impact  loading  situations. 

Finite  Element  Analyses 

Equation  (3-1)  used  in  this  study  in  conjunction  with  the  measured  impact  force  for 
estimating  ILSS  implicitly  assumes  that  the  specimen  is  in  quasi-static  equilibrium 
neglecting  any  dynamic  effects.  In  order  to  verify  this  assumption  Finite  Element 
Analysis  (FEA)  of  the  impact  test  was  performed.  The  analysis  was  simplified  by 
utilizing  the  symmetry  of  the  specimen  and  therefore  only  one  half  of  the  SBS  specimen 
was  modeled  as  shown  in  Figure  3-2. 

In  addition,  the  specimen  was  homogenized  as  an  orthotropic  elastic  material  with 
properties  given  in  Table  3-1  [Xiao  et  al.  2007].  The  use  of  fully  elastic  model  removes 
any  material  rate  dependency  and  focuses  the  study  on  the  effect  of  inertia  and  wave 
mechanics.  The  FEA  study  included  three  rates  of  dynamic  impact  as  well  as  quasi¬ 
static  loading  to  determine  how  loading  rate  affects  the  outcome  of  the  analysis. 

The  commercial  FE  software  codes,  Abaqus  Standard®  and  Abaqus  Explicit® 
were  used  for  the  quasi-static  and  dynamic  loading,  respectively.  The  specimen  was 
modeled  using  four-node  plane  strain  elements  with  reduced  integration  (CPE4R)  and  a 
thickness  set  at  20  mm.  Approximately  32,000  elements  were  used  for  both  static  and 
dynamic  cases.  The  indenter  and  the  support  were  modeled  as  rigid  bodies. 

Dynamic  loading  was  simulated  using  three  different  impact  velocities,  1 1 , 22  and 
33  m/s.  The  maximum  impact  velocity  in  the  experiments  corresponds  to  the  lowest 
rate,  1 1  m/s. 

The  velocity  profiles  for  the  three  impact  simulations  are  depicted  in  Figure  3-3. 
The  shape  of  the  profiles  was  chosen  such  that  they  closely  resemble  the 
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experimentally  measured  profile  discussed  below.  The  duration  of  each  profile  was 
modeled  to  result  in  a  final  displacement  of  2.5  mm.  This  displacement,  which  can  be 
determined  from  the  area  under  the  \/-f  diagram,  is  consistent  with  experimental  results 
[Walter  et  al.  2012].  Both  the  impact  force  (reaction  force  on  the  indenter)  and  the 
support  reaction  were  calculated  for  each  time  increment.  In  the  FEA  the  displacement 
of  the  rigid  indenter  was  controlled  by  the  velocity  profile  up  to  2.5  mm  under  the 
displacement  control.  The  time  step  for  explicit  analysis  was  calculated  automatically  by 
the  FE  program  and  it  was  in  the  order  of  10  ®  seconds. 

Results 

The  impact  force  was  determined  from  the  reaction  on  the  indenter.  The  impact 
force-displacement  curves  for  various  impact  velocities  are  shown  in  Figure  3-4.  The 
deviation  of  the  dynamic  loading  curves  from  the  quasi-static  results  indicates  the 
effects  of  inertia  on  the  load-displacement  relation.  From  these  curves  it  was  observed 
that  up  to  22  m/s,  the  effect  of  inertia  on  the  results  may  be  neglected.  At  33  m/s  this 
effect  becomes  much  more  significant  reducing  the  accuracy  of  the  results.  It  is  clear 
from  these  curves  that  inertia  limits  the  loading  rates  to  approximately  22  m/s. 

Next  the  reaction  force  determined  at  the  support  is  compared  with  the  impact 
force  at  the  indenter  to  verify  quasi-static  equilibrium.  The  results  for  the  each  of  the 
dynamic  loading  simulations  are  shown  in  Figure  3-5  (note  that  the  results  for  quasi 
static  loading  were  not  analyzed  as  it  is  in  static-equilibrium). 

If  the  specimen  is  in  equilibrium  the  reaction  forces  at  the  indenter  and  the  total 
force  at  the  supports  should  be  equal.  From  Figure  3-5A  it  can  be  noted  that  equilibrium 
is  satisfied  during  impact  at  a  velocity  of  1 1  m/s.  At  22  m/s  impact  velocity  (Figure  3-5B), 
the  two  forces  varied  at  the  beginning  of  impact  event  however  after  approximately  0.5 
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mm  of  deflection,  equilibrium  is  reached.  It  is  observed  from  Figure  3-5C  that 
equilibrium  is  never  fully  satisfied  when  loaded  at  a  rate  of  33  m/s. 

Next  the  through  thickness  transverse  shear  stress  profile  was  examined  for  the 
quasi-static,  1 1  m/s  and  22  m/s  results.  Since  the  results  from  33  m/s  impact  were 
previously  shown  to  be  invalid,  due  to  a  state  of  non-equilibrium,  it  was  not  analyzed  for 
stresses.  These  profiles  were  determined  at  a  plane  midway  between  the  support  and 
indenter  (as  shown  in  Figure  3-2).  The  purpose  of  this  comparison  is  to  evaluate  the 
effect  of  loading  rate  on  the  response  of  the  sample  and  compare  the  FEA  profiles  with 
theoretical  profiles.  The  profiles  were  determined  at  a  displacement  of  1 .5  mm  and  are 
presented  in  Figure  3-6. 

The  theoretical  profiles  are  determined  using  the  indenter  force  and  classical 

mechanics  result  for  transverse  shear  stresses  r  =  where  P/2  is  the  shear 

Ih 

force,  Q  is  first  moment  of  area  ,  /  is  moment  of  inertia  and  h  is  the  thickness,)  which  is 
the  basis  for  Eq.  (3-1).  The  FEA  profiles  are  shown  as  solid  lines  while  theoretical 
profiles  are  shown  as  dashed  lines.  It  is  apparent  from  Figure  3-6A  that  there  is  only  a 
slight  deviation  between  quasi-static  and  dynamic  shear  stress  profiles.  This  deviation  is 
shown  much  more  clearly  In  an  enlarged  view  In  Figure  3-6B.  The  theoretical  results 
over-predict  the  FEA  shear  stress  (about  8%  for  quasi-static  and  1 1  m/s  impact  velocity 
case,  and  11%  for  22m/s  impact  velocity  case).  The  shear  stress  results  for  quasl-static 
and  1 1  m/s  Impact  velocity  cases  were  quite  similar;  but  the  results  for  22  m/s  were 
about  3%  higher  than  the  quasi-static  results.  As  mentioned  previously  the  actual  shear 
stress  is  often  much  different  then  what  is  predicted  using  equation  (1)  and  therefore  is 
often  used  to  compare  different  architectures.  These  results  show  that  it  is  possible  to 
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compare  results  from  static  tests  as  well  as  dynamic  tests  at  rates  up  to  1 1  m/s  with  little 
deviation  and  rates  of  22  m/s  with  moderate  error. 

From  the  above  discussion  it  is  clear  that  at  an,  or  below,  impact  velocity  of  1 1  m/s 
the  specimen  may  be  considered  to  be  in  equilibrium,  and  the  use  of  equation  (3-1)  for 
estimating  the  ILSS  of  the  material  could  be  considered  valid. 
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Figure  3-1 .  Short  beam  shear  test. 
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Figure  3-2.  2D  FE  Model  for  the  SBS  test. 
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Figure  3-3.  2D  FE  Model  for  the  SBS  test. 
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(a)  (b) 

Figure  3-5.  Load-displacement  curves  at  the  indenter. 
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Figure  3-6.  Shear  stress  profiles.  A)  Shear  stress  at  the  location  of  midplane  between 
indenter  and  support.  B)  Magnified  view. 


Table  3-1 .  Material  properties  of  the  DSBS  specimen  for  FEA  [Xiao  et  al  2007] 


Material  property 

Value 

Density  (kg/m'*) 

2300 

£i=£2(GPa) 

27.5 

£3  (GPa) 

11.8 

Gi3=G23(GPa) 

2.14 

Gi2(GPa) 

2.9 

Vi3=V23 

0.4 

V12 

0.11 
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CHAPTER  4 

ANALYSIS  OF  DELAMINATION  IN  3D  WOVEN  COMPOSITES  UNDER  QUASI¬ 
STATIC  LOADING  USING  FINITE  ELEMENTS 

Delamination  initiation  and  propagation  in  plain  woven  and  3D  orthogonai  woven 
composite  iaminates  under  three-point  bending  were  anaiyzed  using  finite  element 
analysis.  Both  singie  and  doubie  z-yarn  3D  woven  composites  were  considered.  The 
models  were  guided  by  experimental  observations  of  short  beam  shear  (SBS)  tests  of 
specimens  of  same  materiai  systems.  A  series  of  mechanisms  including  creation  of 
transverse  matrix  cracks  and  their  interaction  with  delamination  were  modeied 
discretely.  The  force-displacement  curves  obtained  from  the  FE  simulation  were 
compared  with  those  from  experiments.  Further  parametric  studies  were  conducted  to 
understand  the  effect  of  z-yarns  in  the  3D  orthogonai  woven  composites  on 
delamination.  The  resuits  from  the  FE  simuiations  reveai  that  z-yarns  in  the  3D  woven 
composites  can  play  a  substantial  role  to  impede  the  propagation  of  interiaminar  cracks, 
thus  can  remarkably  enhance  damage  resistance  and  tolerance  of  the  composite. 

Finite  Eiement  Modei  for  3D  Woven  Composites 

We  selected  three  different  types  of  specimens  -  3D  woven  composites  with  single 
z-yarn  and  double  z-yarn  (SY3W  and  DY3W),  and  2D  plain  woven  laminated  composite 
(2DPL)  -  for  this  study.  The  effect  of  z-yarns  and  inherent  fracture  toughness  of  the 
materiai  on  damage  toierance  was  investigated  by  comparing  the  SY3W  and  DY3W  to 
the  2DPL  specimens  and  varying  inherent  interlaminar  fracture  toughness  of  the  2DPL 
in  a  parametric  study. 

It  is  weil  known  that  matrix  cracks  precede  deiamination  when  a  iaminated 
composite  beam/piate  is  subjected  to  quasi-static  indentation  or  impact  load.  Walter  et 
al.  [Walter  et  al.  2010]  performed  short  beam  shear  tests  on  3D  woven  glass/epoxy 
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composites  to  understand  the  damage  initiation  and  propagation  mechanisms  under 
impact  loads.  The  experiments  showed  that  matrix  cracks  initiated  in  resin  pockets 
between  the  z-yarns  and  also  in  the  weft  yarns  parallel  to  the  y-axis  (Figure  4-1 ). 

Optical  micrographs  also  showed  short  delaminations  and  crushing  failure  beneath 
the  indenter.  However,  these  failures  were  highly  localized,  but  the  delaminations 
emanating  from  matrix  cracks  significantly  changed  the  global  response  of  the 
specimens.  In  fact  the  interaction  of  matrix  cracks  and  delaminations  is  a  common 
feature  typically  observed  in  cross-ply  laminated  composites  comprised  of  0-  and  90- 
layers.  Woven  composites,  which  consist  of  warp  yarns  (0°)  and  weft  yarns  (90°),  are 
similar  to  cross-ply  laminates,  and  understanding  of  damage  mechanism  in  the  latter 
can  be  useful  for  woven  composites.  Although  continuum  damage  models  are  able  to 
capture  the  non-linear  behavior  of  composites  in  certain  situations,  they  are  not  suitable 
for  the  discrete  characteristics  of  damage  process  observed  in  the  present  study 
[Wisnom  2010].  So  the  interaction  between  matrix  cracks  and  delaminations  was 
studied  and  used  for  modeling  the  damage  discretely  [Hallett  et  al.  2008,  Zhou  et  al. 
2010]. 

Cohesive  elements  are  widely  used  for  modeling  of  delaminations  and  matrix 
cracks  because  of  their  versatility  in  dealing  with  damage  initiation  and  evolution  without 
defining  any  pre-crack  [Nishikawa  et  al.  201 1 ,  Okabe  et  al.  2012].  Apart  from 
delamination  which  mostly  occurs  at  the  interface  between  layers,  the  location  and 
direction  of  matrix  crack  may  be  arbitrary.  The  extended  finite  element  method  (X-FEM) 
[Belytschko  et  al.  2009]  or  the  augmented  finite  element  method  (A-FEM)  [Ling  et  al. 
2009]  has  been  used  to  deal  with  these  arbitrary  discontinuities  due  to  cracks  in  the 
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material  system.  However,  it  is  difficult  to  predict  exactly  the  creation  of  multiple  cracks 
in  complex  material  systems  such  as  woven  composites  as  observed  in  the  experiment 
without  considering  uncertainties  in  geometric  configurations,  internal  defects  and 
material  properties.  Since  our  interest  in  the  present  study  is  focused  on  the 
investigation  of  delamination  rather  than  exact  prediction  of  the  intra-laminar  damage 
attributes  in  the  given  material  system,  we  use  the  results  from  experiments  as  a  guide 
to  predict  the  location  of  matrix  cracks.  The  results  from  the  SBS  tests  conducted  by 
Walter  et  al.  [Walter  et  al.  2010]  made  it  possible  to  predict  the  location  of  matrix  cracks 
so  that  cohesive  elements  could  be  used  in  modeling  of  matrix  cracks  as  well  as 
delaminations.  The  potential  locations  for  matrix  cracks  were  assumed  based  on  the 
SBS  tests  and  found  from  principal  stress  directions  for  simplicity,  which  will  be 
described  in  detail  in  the  following  sections. 

As  a  result  of  FE  simulation,  the  force  versus  displacement  at  the  indenter  were 
obtained  for  the  3D  woven  composites  and  used  to  investigate  the  effect  of  z-yarns  as 
well  as  the  interlaminar  fracture  toughness.  The  relationship  between  damage  patterns 
and  the  global  response  was  found  and  will  be  explained  also. 

Geometry  and  Boundary  Conditions 

The  microstructure  of  the  specimens  and  their  boundary  conditions  in  the  FE 
models  were  similar  to  that  in  previous  experiments  [Walter  et  al.  2010].  The  geometric 
parameters  of  the  DY3W  such  as  yarn  dimensions  and  yarn  spacing  were  determined 
from  the  micrographs  of  3D  orthogonal  woven  composite  [Walter  et  al.  2010]  (Table  4- 
1).  As  shown  In  Figures  4-1  and  4-2  the  cross-sections  of  the  yarns  were  treated  as 
rectangles. 
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Taking  advantage  of  symmetry  conditions,  only  a  portion  of  the  specimen  is 
modeled.  Due  to  symmetry  about  the  yz-plane  passing  through  the  mid-span,  only  one- 
half  of  the  specimen  in  the  x-direction  is  modeled.  Since  the  specimen  is  assumed  to  be 
under  plane  strain  normal  to  the  y-axis,  one  unit-cell  in  the  width  direction  (y-direction) 
should  be  sufficient.  Further  reduction  can  be  made  as  the  unit-cell  is  symmetric  about 
its  mid-plane  parallel  to  the  xz-plane.  Thus,  only  one  half  of  the  unit-cell  needs  to  be 
considered  as  shown  in  Figure  4-1 B.  The  indenter  and  support  of  the  SBS  test  setup 
were  modeled  as  rigid  bodies. 

The  span  of  the  beam  (distance  between  the  supports)  in  the  FE  model  was  40 
mm.  The  thickness  of  z-yarn  (tz)  in  the  SY3W  was  0.4  mm  while  other  geometric 
parameters  remained  the  same  as  those  in  the  DY3W  including  the  configuration  of 
crowns  of  z-yarns  in  order  to  investigate  the  effect  of  the  thickness  of  z-yarn,  hence  the 
volume  fraction  of  z-yarns,  on  delamination  of  the  3D  woven  composite.  Since  the 
vertical  part  of  z-yarns  play  a  significant  role  to  interlaminar  crack,  this  approach  might 
be  reasonable.  In  addition,  the  regions  created  by  reducing  the  dimension  of  the  z-yarn 
were  replaced  with  matrix  since  the  micrograph  from  the  experiments  showed  more 
matrix  region  in  the  SY3W  [Walter  et  al.  2010]. 

Plane  strain  boundary  conditions  were  assigned  on  the  front  and  back  surface  of 
3D  woven  composites  (x-z  planes  in  Figure  4-1)  as  well  as  the  symmetry  boundary 
conditions  on  the  symmetry  plane  such  that  Ux=0  (left  y-z  plane  in  Figure  4-1).  Eight- 
node  brick  elements  were  used  in  modeling  the  yarns  and  matrix  phase  and  eight-node 
cohesive  elements  were  used  for  damage  modeling. 
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Damage  Model  and  Modeling  Strategy 

During  the  beginning  stage  of  SBS  tests  3D  woven  composite  specimens  undergo 
elastic  deformation.  As  the  load  increases,  matrix  cracks  and  delaminations  occur  in  the 
specimen. 

The  cohesive  element  associated  with  a  bi-linear  traction-separation  law  was 
selected  for  modeling  of  both  matrix  cracks  and  delamlnatlon.  This  damage  law  enables 
the  traction  between  two  surfaces  to  be  expressed  in  terms  of  relative  displacement  and 
stiffness: 


(7  =  {\-D)Kd 

0  ,d  <dQ 

d ^  ^d  —  d^  ^ 


D  = 


j  I  j  \,dQ<d<d 

d ydf  -d^\ 

1  ,df  <d 


f 


where  cr  is  the  traction,  ^is  stiffness,  Dis  damage  variable,  r/ is  displacement,  i/gisthe 


displacement  at  damage  Initiation  and  df  is  the  final  displacement  [Camanho  et  al. 
2003]. 

Quadratic  stress-based  failure  criterion  and  mixed  mode  strain  energy  release  rate 
criterion  were  used  for  damage  initiation  and  propagation,  respectively,  for  both  the 
transverse  crack  and  delamination: 


where  and  are  interfacial  normal  and  shear  strength,  and  G/c  and  Guc  are  the 


Mode  I  and  Mode  II  critical  energy  release  rates.  Since  the  locations  of  cohesive 
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elements,  which  would  act  as  potential  crack  paths,  should  be  defined  a  priori,  it  is 
necessary  to  predict  accurately  the  locations  of  expected  cracks.  The  prediction  for 
delamination  is  straightforward  as  it  occurs  between  adjacent  plies  or  layers,  while  that 
for  matrix  crack  is  complex.  Based  on  experimental  observations  [Walter  et  al.  2010] 
several  assumptions  have  been  made  to  determine  the  locations  of  potential  cracks  for 
the  efficient  use  of  cohesive  elements.  The  assumptions  are  as  follows: 

i)  Matrix  cracks  occur  both  within  the  fill  tows  and  in  the  matrix  pocket  at  the 
tensile  bottom  layer. 

ii)  Only  a  single  matrix  crack  is  created  and  allowed  to  evolve  along  the  principal 
stress  direction. 

iii)  The  crack  path  is  a  straight-line.  The  above  assumptions  are  consistent  with 
experimental  observations. 

iv)  The  effect  of  a  matrix  crack  at  the  compressive  top-most  layer  on  delamination 
can  be  negligible. 

Formation  of  matrix  cracks  on  the  top  layers  of  the  beam  needs  special  treatment. 
In  the  beginning  of  loading  the  top  side  is  under  compression  and  it  is  observed  no 
matrix  cracks  form.  However,  after  delaminations  initiate,  there  is  a  redistribution  of 
stresses  and  the  delaminated  top  layers  develop  tensile  stresses  and  they  become  sites 
for  initiation  of  matrix  cracks.  These  locations  are  identified  from  a  preliminary  FE 
analysis.  To  begin  with,  the  principal  directions  at  the  centers  of  fill  tows  in  the  tensile 
region  using  the  specimen  without  any  damage  (Figure  4-3A)  were  found  and  used  for 
further  simulation  to  seek  the  principal  stress  direction  over  the  entire  region  as  seen  in 
Figure  4-3B.  In  the  following  FE  model,  other  cohesive  elements  for  delaminations  were 
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placed  between  layers  with  the  cohesive  element  for  matrix  cracks  in  the  tensile  region 
placed  along  the  principal  direction  found  in  the  previous  step,  thus  the  principal 
directions  in  the  compressive  region  after  delamination  could  be  found.  In  the  second 
FE  analysis  using  the  specimen  of  Figure  4-3B  the  load  was  applied  until  the  moment 
that  delamination  occurred  so  that  redistributed  stress  filed  could  be  obtained.  In  order 
to  determine  the  principal  stress  directions  in  the  weft  yarns,  stress  values  over  the  all 
elements  of  the  weft  yarns  were  checked  and  appropriate  principal  stress  directions 
incurring  tensile  failure  were  chosen.  The  principal  directions  determined  are  shown  in 
Figure  4-3.  Note  that  the  angles  shown  in  Figure  4-3  are  measured  from  the  z-axis.  As  a 
result  of  consecutive  finite  element  simulations,  the  final  configuration  of  Figure  4-3C 
could  be  obtained  and  became  the  final  FE  model  for  the  study  of  delamination  damage 
of  the  3D  woven  composite  materials. 

Material  Property  and  Data  Selection 

Material  properties  of  the  3D  woven  composites  used  in  the  analyses  are  listed  in 
Table  4-2.  Elastic  properties  of  the  yarns  were  computed  using  the  micromechanics 
formulas  developed  by  Chamis  et  al.  [Chamis  et  al.  2007].  The  material  system  was  S2 
glass  fiber/  SC-15  epoxy  with  62%  fiber  volume  fraction.  The  strength  and  critical 
energy  release  rate  values  for  the  cohesive  elements  were  selected  from  the  literature. 
The  properties  of  cohesive  elements  for  matrix  crack  and  delamination  were  assumed 
to  be  same. 

As  Mode  II  shear  strength  value  for  the  cohesive  elements,  the  apparent  ILSS 
value  reported  by  Walter  et  al.  [Walter  et  al.  2010]  was  chosen  (Table  4-2).  This  value 
was  much  lower  than  the  range,  55  to103  MPa,  given  by  AGY  manufacturer  of  S2  glass 
yarn,  [AGY].  The  ILSS  could  be  much  lower  due  to  the  voids  and  the  type  of  binder 
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used  in  the  manufacturing  process  [Wisnom  et  al.  1996,  Tanoglu  et  al.  2001].  The 
specimens  test  by  Walter  et  al.  were  manufactured  using  vacuum  assisted  resin  transfer 
molding  (VARTM)  process  [Walter  et  al.  2010],  the  ILSS  value  could  be  lower  because 
of  inherent  defects  such  as  voids  involved  during  VARTM  process.  Based  on  the  Mode 
II  shear  strength  determined  for  the  cohesive  element,  a  Mode  I  normal  strength  value 
was  estimated  appropriately  as  listed  in  Table  4-2 

In  addition,  the  critical  energy  value  for  the  cohesive  element  could  be  determined 
using  available  sources  such  as  a  Mode  II  fracture  toughness  value  of  quasi-isotropic 
and  cross-ply  composite  with  SC-15  reported  by  Huang  [Huang  2008],  which  was  about 
2000  N/m,  as  well  as  a  Mode  I  fracture  toughness  of  SC-15  epoxy  provided  by  Applied 
Poleramic,  the  provider  of  SC-15  epoxy,  which  was  1000  N/m  [Applied  Poleramic].  A  bi¬ 
linear  traction-separation  law  controlling  the  behavior  of  the  cohesive  element  can  be 
depicted  by  a  triangle  as  shown  in  Figure  4-4. 

If  the  critical  displacement  denoted  by  ‘C’  in  Figure  4-4  remains  the  same,  the 
critical  energy  will  be  changed  (from  the  area  of  ‘OBD’  to  the  area  of  ‘OAD’)  with  the 
change  of  the  strength  (from  ‘B’  to  ‘A’).  This  approach  can  explain  the  fact  that  fracture 
toughness  tends  to  decrease  due  to  the  lowered  strength  [Madhukar  et  al.  1992]. 
Otherwise  the  lowered  strength  value  can  enhance  damage  resistance  by  increasing 
the  critical  distance  of  cohesive  element  under  the  assumption  of  same  fracture 
toughness  value  (from  the  area  of  ‘OBD’  to  the  area  of  ‘OAE’)  and  might  not  be  practical 
for  this  study.  With  these  considerations,  the  critical  energy  value  for  the  cohesive 
element  was  determined  as  listed  in  Table  4-2. 
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Finite  Element  Model  for  2D  Plain  Woven  Laminates 

A  plain  woven  laminate  (2DPL)  with  stacking  sequence  (o^'^/45^'^)2/o^'^  5'^^® 

analyzed.  The  2DPL  was  chosen  for  comparison  with  3D  woven  composites  thus 
investigating  the  effect  of  z-yarns  on  delamination  behavior.  For  the  sake  of  simplicity  of 
the  FE  simulation,  each  plain  woven  lamina  was  homogenized  as  an  orthotropic 
material  and  one-half  of  the  beam  was  modeled  using  eight-node  plane  strain  elements 
(Figure  4-5).  The  material  properties  of  the  homogenized  plain-woven  layer  are  listed  in 
Table  4-3  [Xiao  et  al.  2007]. 

The  material  properties  for  ±45^'^  could  be  obtained  using  coordinate 
transformation. 

The  creation  of  matrix  cracks  and  placement  of  cohesive  elements  in  2DPL 
specimens  were  similar  to  the  procedures  used  for  3D  woven  laminates.  Our  approach 
is  confirmed  by  the  experimental  observation  that  a  single  matrix  crack  tended  to  occur 
inside  weft  yarns  of  the  O/OO^"^  [Walter  et  al.  2010].  Flowever,  the  tensile  cracks  in  the 
matrix  pocket  of  outer  most  layer  seen  in  the  3D  woven  composites  were  not  observed 
in  plain  woven  laminates  and  hence  they  were  not  introduced  in  2DPL  specimens. 
Furthermore,  transverse  shear  stress  vanishes  at  the  free  surface.  Therefore  cohesive 
elements  were  merely  placed  in  the  middle  regions  of  the  specimen  only.  Cohesive 
elements  for  possible  matrix  cracks  were  placed  in  the  three  inner  homogenized  layers 
of  0/90^'^as  shown  in  Figure  4-6.  Principal  stress  directions  were  found  in  locations 
corresponding  to  the  centers  of  fill  tows  through  same  procedure  as  3D  woven 
composites  considering  tensile  region  and  compressive  region  after  delamination 
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individually.  The  possible  matrix  crack  directions  used  in  this  simulation  are  shown  in 
Figure  4-6. 

A  parametric  study  varying  Mode  I/Mode  II  interlaminar  fracture  toughness 
(370/830  N/m,  1000/2000  N/m  and  1000/3300  N/m)  was  performed  in  order  to 
investigate  the  effect  of  interlaminar  fracture  toughness  on  the  global  behavior  and 
damage  pattern  of  the  2DPL  specimens.  These  studies  are  referred  to  as  Case  I,  Case 
II  and  Case  III  (Table  4-4),  respectively. 

Results  and  Discussions 
3D  Woven  Composites  (SY3W  and  DY3W) 

The  force  and  displacement  at  the  indenter  were  recorded  during  the  FE 
simulations.  The  force  at  a  given  displacement  can  be  a  proper  measure  to  examine  the 
macroscopic  behavior  the  specimen.  The  results  for  both  SY3W  and  DY3W  are  shown 
in  Figure  4-7. 

Two  force-displacement  curves  from  the  tests  [Walter  et  al.]  are  also  plotted  for 
comparison.  Note  the  two  experimental  curves  of  the  SY3W  type  show  slightly  different 
behavior  after  damage  initiation.  The  solid  lines  indicate  the  result  from  the  FEM.  Since 
the  same  properties  for  the  cohesive  element  were  used  and  most  configurations 
remained  the  same  except  for  the  thickness  of  z-yarns  in  the  simulation  it  can  be 
inferred  that  the  difference  between  two  FEM  simulations  (solid  curves  in  Figure  4-7) 
resulted  from  the  effect  of  z-yarns  on  delamination.  The  peak  load  value  of  DY3W  is 
higher  than  that  of  SY3W.  This  reveals  that  z-yarns  can  delay  the  damage  initiation 
point.  Besides  the  forces  beyond  the  peak  loads  between  two  cases  are  quite  different. 
The  SY3W  showed  that  an  abrupt  force-drop  right  after  the  peak  loads  whereas  the 
force-drop  in  DY3W  is  not  noticeable  or  fairly  small.  Additionally  the  SY3W  has  another 
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peak  load  followed  by  sudden  force  drop  again  which  does  not  occur  in  the  DY3W. 
These  characteristics  agree  well  with  experimental  results  shown  in  Figure  4-7. 

The  variation  of  strain  energy  and  damage  energy  as  a  function  of  indenter 
displacement  obtained  from  the  FE  analyses  are  shown  in  Figure  4-8. 

The  strain  energy  represents  the  recoverable  elastic  energy  stored  in  the 
specimen  at  any  given  instant.  The  damage  energy  is  the  unrecoverable  energy 
dissipated  by  the  cohesive  elements.  It  is  seen  that  the  sudden  force-drop  in  the  load- 
deflection  diagram  corresponds  to  sudden  increase  in  damage  energy.  One  can  note 
that  the  damage  energy  in  the  DY3W  specimen  increases  steadily  without  any  abrupt 
damage  process  or  loss  of  stiffness.  From  these  figures  it  is  clear  that  z-yarn  can 
enhance  damage  resistance  and  tolerance.  However,  optimum  amount  of  z  yarn  should 
be  used  so  that  there  is  no  deterioration  in  the  composite  properties  [Rao  et  al.  2009  ]. 
Plain  Woven  Laminated  Composites  (2DPL) 

Unlike  3D  woven  composites,  only  one  configuration,  but  with  various  cohesive 
parameters,  was  considered  in  studying  of  the  2DPL.  Although  the  same  properties  for 
cohesive  element  as  3D  woven  composites  were  Initially  used  for  the  2DPL, 
interlaminar  fracture  toughness  values  seems  to  be  different  since  woven  laminated 
composites  have  higher  fracture  toughness  than  unidirectional  fiber  composite  due  to 
undulation  of  yarn  or  crimp  [Kalarikkal  et  al.  2006].  The  force-displacement  curves  of 
both  SY3W  and  2DPL  obtained  from  the  tests  [Walter  et  al.  2010]  showed  different 
behavior  during  unloading.  The  force  drop  of  the  SY3W  was  very  rapid  while  that  of  the 
2DPL  was  gradual.  This  could  be  explained  by  the  higher  fracture  toughness  values  of 
woven  laminates  compared  to  tape  laminates.  Out  of  the  three  sets  of  interlaminar 
fracture  toughness  values  used  the  force-displacement  curve  of  Case  III  with  the 


69 


highest  fracture  toughness  exhibited  load-displacement  behavior  similar  to  experimental 
results  as  shown  in  Figure  4-9. 

Lower  fracture  toughness  values  (cases  I  and  II)  resulted  in  different  damage 
pattern  and  produced  different  load-displacement  curves. 

The  force  versus  displacement  curves  at  the  indenter  with  respect  to  three 
different  interlaminar  fracture  toughness  values  are  plotted  in  Figure  4-9.  The  difference 
among  the  force-displacement  curves  due  to  different  fracture  toughness  values  is 
apparent.  As  the  fracture  toughness  values  were  increased,  the  peak  loads  were 
increased  and  the  rate  of  force  drop  just  beyond  the  peak  load  became  moderate.  If  the 
interlaminar  cracks  propagate  mainly  without  other  delaminations  the  forces  will 
decrease  rapidly.  As  seen  in  Figure  4-9,  the  rapid  force  drops  in  the  Case  I  and  Case  II 
were  attributed  to  the  damages  that  occurred  with  very  small  increment  of  displacement 
and  can  be  identified  with  the  damage  energy  variation  in  Figure  4-10.  Sudden  increase 
of  damage  energy  was  accompanied  by  sudden  loss  of  strain  energy  as  well.  But  the 
damage  energy  of  Case  III  increases  gradually  along  with  the  increase  of  the  strain 
energy. 

These  variations  of  global  responses  are  highly  dependent  on  how  damage 
evolves.  By  observing  the  history  of  damage  development,  it  is  possible  to  gain  insight 
into  understanding  delamination  damage  behavior  of  composites.  Once  delamination 
occurs  from  the  matrix  crack  it  will  propagate  and  reach  the  free  end  of  specimen  (b-1 
and  b-2  in  Figure  4-11). 

The  propagation  of  the  crack  to  the  free  end  causes  the  force  to  decrease 
abruptly.  After  a  complete  delamination  developed  over  the  region  from  matrix  crack  to 
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the  end  of  the  specimen,  the  force  increases  with  displacement  and  strain  energy  is 
stored  until  the  state  of  stress  in  another  interlaminar  region  satisfies  the  criterion  for 
crack  initiation  and  propagation.  Then  the  force  will  decrease  suddenly  again  as  shown 
in  the  Case  I  of  Figures  4-9,  4-10  and  Figure  4-11 C.  Case  II  showed  that  the  resistance 
to  crack  propagation  increased  with  the  interlaminar  fracture  toughness.  Further 
increase  in  interlaminar  fracture  toughness  (Case  III)  totally  changed  the  tendency  of 
crack  propagation. 

High  resistance  against  cracks  kept  the  interlaminar  delamination  from 
propagating.  Before  a  delamination  reached  the  end  of  specimen,  another  delamination 
occurred.  This  caused  a  gradual  force  drop  instead  of  an  abrupt  drop  iqn  Cases  I  and  II 
(Figures  4-9,  4-10  and  4-13).  The  relationships  between  the  damage  pattern  in  Figures 
4-10  through  4-12  and  the  force-displacement  curves  in  Figure  4-9  are  similar  to  the  test 
results  of  Walter  [Walter  201 1  ]. 

It  could  be  observed  that  the  crack  propagation  of  the  SY3W  reached  the  end  of 
the  specimen  and  resulted  in  the  force-displacement  curve  in  Figure  4-7,  while  2DPL 
and  DY3W  didn’t  allow  delamination  to  propagate  to  the  end  of  the  specimen  and  other 
delaminations  occurred  at  the  same  time  they  resulted  in  gradual  force  drop  as  shown 
in  Figure  4-14.  It  should  be  noted  that  Figure  4-14  has  been  obtained  from  a  dynamic 
test  but  exhibited  similar  damage  pattern  as  quasi-static  test  demonstrating  that  quasi¬ 
static  indentation  tests  can  provide  insight  into  damage  development  due  to  impact 
loading  in  laminated  3D  composites. 

From  these  results,  it  is  possible  to  identify  that  the  role  of  z-yarns  on  interlaminar 
delamination.  Z-yarns  provide  a  constraint  to  reduce  the  relative  displacement  between 
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two  layers  thus  enhancing  the  resistance  against  crack  propagation.  Thus  the  z-yarns 
result  in  higher  apparent  interlaminar  fracture  toughness. 

Summary  and  Conclusion 

The  effect  of  z-yarns  on  damage  tolerance  of  3D  woven  composites  was  studied 
using  the  SBS  tests.  It  is  found  that  tensile  or  transverse  matrix  cracks  develop  in  the 
matrix  pocket  between  z-yarn  crowns  and  the  center  of  wept  yarns.  Cohesive  elements 
associated  with  bi-linear  damage  model  were  placed  along  the  principal  stress  direction 
for  the  matrix  cracks.  These  elements  are  used  in  modeling  of  interlaminar  delamination 
as  well.  The  RVE  model  and  simplified  plain  strain  FE  model  provided  good  results  to 
understand  the  role  of  z-yarns  and  inherent  interlaminar  fracture  toughness. 

From  the  FE  simulation  of  SY3W  and  DY3W  specimens,  it  is  found  that  proper 
amount  of  z-yarns  can  enhance  damage  resistance  and  tolerance.  This  conclusion 
agrees  well  with  the  result  from  the  end  notch  flexure  (ENF)  specimen  of  3D  woven 
composites  by  Pankow  et  al.  [Pankow  et  al.  2011].  It  was  found  that  the  z-yarns  prevent 
two  neighboring  layers  from  being  separated.  In  the  case  of  2DPL  crimp  or  undulation  of 
the  yarns  seems  to  interrupt  the  propagation  delamination  crack  providing  better 
damage  tolerance  than  single  z-yarn  3D  woven  composites.  Flowever,  double  yarn  3D 
woven  composites  (DY3W)  have  superior  damage  tolerance  compared  to  2DPL  and 
SY3W.  As  it  turned  out,  the  prevention  of  crack  growth  along  an  entire  interface  is  a  key 
element  to  enhance  damage  tolerance.  When  damage  occurs  in  different  layers  before 
the  delamination  propagates  in  the  entire  specimen,  more  energy  is  dissipated  in 
various  interlayer  damage  mechanisms  increasing  the  damage  tolerance  of  the 
laminated  composite  structure. 
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A  B 

Figure  4-1 . 3D  woven  composites  for  FE  model.  A)  Orthogonally  woven  3D  composite 
with  SBS  setup.  B)  Representative  Volume  Element  (RVE)  for  FE  model. 


A  B 

Figure  4-2.  Geometric  parameter  of  RVE.  A)  left  side  view.  B)  front  side  view. 
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Figure  4-3.  Prediction  of  matrix  cracks.  The  iines  indicate  potential  crack  paths.  A) 
Cross-section  for  exploring  principai  stress  directions.  B)  Principal  stress 
directions  in  the  tensiie  region.  C)  Principal  stress  directions  in  the  entire 
regions. 


74 


Figure  4-4.  Traction-separation  law  of  cohesive  element. 


Figure  4-5.  2D  plain  woven  composite.  A)  Cross-section  of  2D  plain  laminated 
composite.  B)  Homogenized  FE  model. 
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Figure  4-6.  Implementation  of  the  cohesive  element  on  homogenized  2D  plain  woven 
composite. 
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Energy  (J)  Cq'  Force  (kN) 


ure  4-7.  Force-displacement  curves  of  the  single/double  z-yarn  3D  woven 
composites. 


Figure  4-8.  Strain  energy  and  damaged  energy  versus  displacement  curves  of  the 
SY3W  and  DY3W.  SE  and  DE  indicate  strain  energy  and  damage  energy, 
respectively. 
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Figure  4-9.  Force-displacement  curves  of  the  single  and  double  z-yarn  3D  woven 
composites. 


Figure  4-10.  Strain  energy  and  damage  energy  versus  displacement  curves  of  the 

2DPW.  SE  and  DE  indicate  strain  energy  and  damage  energy,  respectively. 
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Figure  4-1 1 .  Damage  evolution  of  baseline  composite-Case  I.  A)  d=0.8mm.  B)  d=1 .1  mm. 
C)  d=1.7  mm. 
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A 


B 


Figure  4-12.  Damage  evolution  of  baseline  composite-Case  II.  A)  d=1mm.  B)  d=1.5mm. 
C)  d=1 .75  mm. 
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Figure  4-13.  Damage  evolution  of  baseline  composite-Case  III.  A)  d=1mm.  B)  d=1.5mm. 
C)  d=1.75  mm. 


80 


ABC 


Figure  4-14.  Delamination  patterns  in  the  composites.  A)  2DPL.  B)  SY3W.  C)  DY3W. 
[Waiter  2012]. 
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Table  4-1 .  Dimensions  of  geometric  model  of  RVE.  The  subscripts  ‘x’,  ’y’  and  ’z’  indicate 
_ warp  yarn,  weft  yarn  and  z  yarn,  respectively,  (unit:  mm) _ 

U  Wz  tx  ty  tz  Wy  tg 


1.8  0.4  0.5  0.6  0.8  2.8  0.075 


Table  4-2.  Material  properties  of  the  constituents  the  3D  woven  composites 

Yarns,  Epoxy  and  Cohesive  element 

Yams  (Transversely  isotropic) 


El  E2  E3  Gi2  G23  G31 


(GPa) 

(GPa) 

(GPa) 

V12 

Vis 

^23  (GPa)  (GPa) 

(GPa) 

54.9 

11.4 

11.4 

0.26 

0.26 

0.29  4.4  4.4 

4.4 

Epoxy 

Cohesive  element 

E 

(GPa) 

2.7 

V 

0.3 

Mode  I/ll 
Strength 
(MPa) 

Mode  I/ll 

23/33 

toughness 

(N/m) 

370/830 

Table  4-3.  Elastic  properties  for  plain  woven  composites 


Ei=E2 

(GPa) 

E3 

(GPa) 

G13-G23 

(GPa) 

Gi2 

(GPa) 

V13-V23 

V12 

27.5 

II.8 

2.14 

2.9 

0.4 

O.II 

Table  4-4. 

Parameters  for  cohesive  element  in  baseline  composite  model 

^max 

Gic 

^max 

Giic 

(MPa) 

(N/m) 

(MPa) 

(N/m) 

Case  I 

23 

370 

33 

830 

Case  II 

23 

1,000 

33 

2,000 

Case  III 

23 

1,000 

33 

3,300 
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CHAPTER  5 

DAMAGE  RESISTANCE  OF  HYBRID  COMPOSITES  WITH  FUNTIONALLY  GRADED 

MATERIALS 

Along  with  the  attempts  to  enhance  material  properties  such  as  strength  and 
fracture  toughness  a  way  to  reduce  stress  or  strain  energy  reiease  rate  (SERR) 
subjected  by  externai  loading  can  be  a  significant  consideration  in  the  design  and 
operation  of  composite  structures.  That  is  because  the  stress  or  strain  energy  release 
rate  values  in  the  material  should  be  less  than  the  strength  or  fracture  toughness  of  the 
material  in  order  to  prevent  failure  of  a  material.  For  laminated  structures,  interiaminar 
shear  stress  can  be  a  proper  measure  to  evaluate  the  damage  resistance  of  materiais. 
As  discussed  in  previous  chapters  the  beam-iike  specimens  subjected  to  quasi-static  or 
low  velocity  impact  ioading  exhibit  a  typical  shear  stress  distribution  wherein  the 
maximum  value  occurs  beneath  the  contact  surface.  Although  the  shear  stresses  at  this 
location  create  deiamination,  they  do  not  propagate  due  to  rapidly  diminishing  shear 
stresses  away  from  the  loading  point.  The  shear  stresses  are  almost  negligible  near  the 
top  surface  of  the  beam  away  from  the  point  of  contact.  There  is  another  location  where 
the  shear  stresses  are  high  enough  to  cause  deiaminations.  This  is  at  the  mid-piane  of 
the  beam  where  the  paraboiic  shear  stress  profiie  attains  a  maximum.  In  fact  this  region 
of  maximum  shear  stress  is  substantiaiiy  iarge.  Hence,  delaminations  can  propagate  all 
along  the  length  of  the  specimen  until  they  reach  the  ends  as  the  shear  force  is  constant 
in  SBS  specimens  extending  the  region  of  maximum  shear  stress.  In  general,  the 
delamination  due  to  an  impact  ioading  initiates  when  interlaminar  shear  stress  exceeds 
interiaminar  shear  strength. 

As  a  method  to  reduce  the  maximum  interiaminar  shear  stress,  hybrid  composites 
comprised  of  more  than  two  materials  or  functionally  graded  materials  (FGM)  with 
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continuously  varying  material  property  can  be  considered.  In  the  present  study  impact 
damage  resistance  of  hybrid  laminated  composites  and  functionally  graded  laminated 
composited  are  investigated.  The  motivation  for  this  approach  comes  from  some  earlier 
studies  on  functionally  graded  beams  [Sankar  2001]  and  hybrid  laminates  [Sankar 
1989].  It  was  found  that  the  transverse  shear  stresses  get  redistributed  in  FG  beams 
and  hybrid  laminates. 

The  specimens  of  hybrid  composites  studied  have  various  combinations  based  on 
the  volume  fraction  of  carbon  fiber  composite  (CFG)  and  glass  fiber  composite  (GFC), 
stacking  sequence  while  the  functionally  graded  composite  have  three  different  aspects. 
The  descriptions  of  those  specimens  are  shown  in  Table  5-1 . 

Effect  of  Hybrid  Composites  on  Reducing  the  Maximum  Shear  Stress 
Finite  Eiement  Modei  for  Maximum  Shear  Stress 

The  shear  stress  profiles  including  maximum  values  in  various  specimens  under 
an  assumed  force  were  found  using  FEM  since  the  evaluation  of  damage  resistance  of 
various  materials  could  be  achieved  by  comparing  the  maximum  interlaminar  shear 
stress.  The  force  in  the  FE  model  was  assumed  as  a  Hertzian  contact  force  given  by 

I  fxY 

p  =  100,1-  —  N/rnm,  |x|  <  3rnm 

V  v3y 

The  FE  model  to  determine  the  shear  stress  profiles  is  shown  in  Figure  5-1 .  Due  to 
symmetry,  one  half  of  the  specimen  was  modeled  using  plane  strain  elements.  The 
dimensions  of  the  half-specimen  was  100  mm  x15.4  mm.  Symmetric  conditions  were 
implemented  at  x=0,  and  the  vertical  displacements  at  x=100  mm  were  constrained. 

The  shear  stress  profiles  at  the  location  of  x=50  mm  (50  mm  from  the  symmetry  plane) 
were  obtained  from  the  FE  model. 
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The  material  properties  used  in  the  FE  model  are  given  in  Table  5-2.  The 
properties  for  the  FGM  were  found  using  the  weighted  average  based  on  volume 
fraction  of  the  CFG  and  the  GFC.  This  approach  was  based  on  the  resent  results 
obtained  by  Banerjee  [Banerjee  2012]  using  micromechanics  methods  for  hybrid 
composites. 

Results  and  Discussion 

Hybrid  composites  with  loading  on  the  harder  (stiffer)  side.  Figure  5-2  shows 
the  shear  stress  profiles  of  hybrid  composites  with  carbon  fiber  laminates  placed  on  the 
loading  side.  In  the  same  figure  shear  stress  results  for  FGM  with  linearly  decreasing 
stiffness  from  the  loading  side  to  the  free  side  are  also  presented.  The  maximum  shear 
stress  increases  with  the  amount  of  carbon  fiber  layer.  The  maximum  shear  stress  value 
of  the  LS36Hybrid  is  almost  same  as  that  of  homogeneous  composite  such  as  the  GFC 
or  the  CFG.  It  should  be  noted  that  the  shear  stress  profile  of  the  GFC  is  exactly  the 
same  as  that  of  the  CFG.  The  maximum  shear  stress  in  the  LSFGM  is  higher  than  that 
of  the  GFC  or  CFG.  Interestingly,  the  location  of  the  maximum  share  stress  is  shifted 
from  the  midplane  toward  the  loading  side,  and  thus  can  affect  the  point  where 
delamination  will  occur. 

Hybrid  composites  with  loading  on  the  softer  side.  Figure  5-3  shows  the  shear 
stress  profiles  of  the  hybrid  composites  with  carbon  fiber  laminates  placed  on  the 
bottom  (free)  side.  As  the  volume  fraction  of  carbon  fiber  is  increased,  the  maximum 
shear  stress  become  higher.  The  FSFGM  has  higher  maximum  shear  stress  than  the 
homogenous  material.  The  location  of  the  maximum  shear  stress  is  shifted  from  the 
middle  toward  the  free  side  direction. 
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Hybrid  composites  stiffened  on  the  top  and  bottom  sides.  Figure  5-4  shows 
the  shear  stress  profiles  of  the  hybrid  composites  with  carbon  fiber  laminates  placed  on 
the  both  loading  and  free  sides.  The  locations  of  the  maximum  shear  stress  values  in  all 
the  hybrid  composites  stiffened  on  both  sides  are  at  the  middle.  The  maximum  shear 
stress  values  of  both  the  hybrid  composites  and  the  FGM  are  lower  than  those  of  the 
composites  made  of  a  single  material.  The  DS18Hybrid  composite  has  the  lowest  value 
while  the  DSFGM  has  the  value  between  the  hybrid  composites  and  the  composites  of 
a  single  material. 

Functionally  Graded  Composites.  Figure  5-5  shows  the  shear  stress  profiles  of 
the  FGMs.  The  DSFGM  has  lower  maximum  shear  stress  level  compared  to 
homogeneous  composites  while  the  LSFGM  and  the  FSFGM  have  higher  maximum 
shear  stress  values.  Like  the  hybrid  composites,  the  locations  of  the  maximum  shear 
stress  in  the  LSFGM  and  the  FSFGM  are  away  from  the  midplane. 

Comparison  of  maximum  shear  stresses.  Figure  5-6  shows  the  comparison  of 
the  maximum  shear  stress  values.  The  hybrid  composite  DS18Hybrid  has  the  lowest 
value  which  is  about  1 1  percent  less  than  the  homogeneous  beam  (GFC  or  the  CFG) 
while  both  the  LSFGM  and  the  FSFGM  have  the  highest  value  which  is  about  2.8 
percent  more.  It  can  be  expected  that  the  specimen  with  lower  shear  stress  level  has 
higher  damage  resistance  because  an  extra  force  or  energy  would  be  required  to  reach 
the  shear  strength  value. 

Delamination  Behavior  of  Hybrid  Composites  under  a  Quasi-static  Indentation 
Finite  Element  Model  for  Delamination 

It  can  be  expected  that  the  hybrid  composite  stiffened  on  both  sides  would  show 
higher  peak  load  over  the  one-material  composite  specimens  since  lower  shear  stress 
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has  more  margin  to  reach  the  strength  value.  Two  FE  models  to  determine  peak  loads 
of  various  specimens  described  in  Table  5-1  under  3  points  bending  were  considered. 
(Figure  5-7C)  One  has  only  delaminations.  (Figure  5-7A)  The  other  has  matrix  cracks 
and  delaminations  (Figure  5-7B).  All  specimens  in  the  FE  models  consist  of  22  layers. 
Cohesive  elements  were  placed  between  layers  as  mentioned  in  Chapter  4  thus 
delaminations  were  allowed  to  be  included  in  the  FE  models.  In  the  case  of  the 
specimen  with  matrix  cracks,  both  the  matrix  cracks  and  the  delaminations  are  modeled 
using  the  cohesive  elements.  The  matrix  cracks  were  placed  at  the  location  of  x=37.5 
mm  from  the  symmetry  line.  Plane  strain  elements  were  used  for  the  half  model  due  to 
symmetry.  Both  the  indenter  and  the  support  were  modeled  using  rigid  bodies.  The 
displacement  of  indenter  was  defined  as  a  boundary  condition  and  the  load 
corresponding  to  each  displacement  was  found.  The  peak  loads  from  both  FE  models 
were  compared  so  as  to  analyze  the  effect  of  matrix  cracks  on  damage  resistance.  It  is 
found  that  the  matrix  cracks  altered  the  layer  in  which  delaminations  occurred. 

Results 

Figure  5-8  shows  the  peak  loads  of  all  hybrid  specimens  considered  in  this  study. 
Figures  5-9  through  5-20  show  the  load-deflection  curves  and  the  damage  patterns  in 
various  specimens.  Obviously  the  specimens  with  matrix  cracks  exhibit  lower  peak 
loads.  This  is  because  there  is  stress  concentration  at  the  tip  of  the  matrix  cracks 
enabling  the  delaminations  to  propagate  at  lower  loads.  As  the  matrix  cracks  were 
introduced  the  location  of  delamination  shifts  from  the  maximum  interlaminar  stress 
location  to  the  interface  that  contains  the  crack  tip.  The  DSISHybrid  composite  has  the 
highest  peak  load  and  hence  highest  damage  resistance.  Detailed  discussions  about 
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the  location  of  delamination  as  well  as  load-displacement  curve  for  individual  specimens 
are  given  in  the  following. 

Figure  5-9  and  Figure  5-10  show  the  load-displacement  curve  and  the 
deformations  after  delamination  of  the  GFC  and  the  CFG,  respectively.  The  location  of 
delamination  changed  from  the  11 interface  (the  middle  of  the  specimen)  to  the  8**^ 
interface.  As  long  as  the  specimen  has  a  symmetric  configuration  delamination  occurs 
at  the  middle  since  the  shear  stress  is  maximum  at  the  center.  However  the  matrix 
cracks  cause  the  stress  field  to  be  redistributed  and  change  the  delamination  site. 

Figures  5-1 1  through  Figure  5-14  show  the  load-displacement  curve  and  the 
deformations  after  delamination  of  the  LS9Hybrid  composite,  the  LSISHybrid 
composite,  the  LS36Hybrid  composite  and  the  LSFGM  composite  (Loading  side  was 
stiffened).  Figure  5-15  through  5-18  show  the  load-displacement  curves  and  the 
deformations  after  delamination  of  the  FS9Hybrid  composite,  the  FS18Hybrid 
composite,  the  FS36Hybrid  composite  and  the  FSFGM  composite  (Free  side  was 
stiffened).  Figure  5-19  shows  the  load-displacement  curves  and  the  deformations  after 
delamination  of  the  DS9Hybrid  and  Figure  5-20  shows  the  load-displacement  curves  of 
other  cases  in  which  both  sides  were  stiffened.  In  the  case  of  the  specimens  without 
matrix  cracks  the  location  of  delamination  can  be  determined  as  the  location  of  the 
maximum  shear  stress.  Therefore  the  location  of  delamination  varied  based  on  the 
configurations.  However,  in  all  cases  with  matrix  cracks  the  delamination  occurred  at 
the  10**^  interface  for  the  case  of  loading  on  the  stiffer  side,  at  the  6^'^  interface  for  the 
case  loading  on  the  softer  side,  and  at  the  8**^  interface  when  both  sides  were  stiffened 
(symmetric  beam). 
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Conclusions 


Hybrid  composites  including  FGM  composites  can  reduce  the  maximum  shear 
stress,  and  thus  can  enhance  damage  resistance.  However  the  un-symmetric  FGM 
composites  results  in  higher  maximum  shear  stress.  The  effects  of  reducing  the 
maximum  shear  stress  become  noticeabie  when  stiff  materials  are  symmetrically  used 
in  both  sides.  The  locations  of  delamination  of  the  doubly  stiffened  specimens  are 
identical  with  the  homogeneous  composite.  However,  the  delaminations  of  the  singly 
stiffened  specimen  occur  at  the  interface  toward  the  stiffened  side.  Matrix  cracks  alter 
the  locations  of  delamination  and  decrease  damage  resistance. 
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Figure  5-1 .  FE  model  to  determine  the  shear  stress  profile 
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Figure  5-2.  Shear  stress  profiles  of  the  hybrid  composites  stiffened  on  the  loading  side 
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Figure  5-3.  Shear  stress  profiles  of  the  hybrid  composites  stiffened  on  the  free  side 
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Figure  5-4.  Shear  stress  profiles  of  the  hybrid  composites  stiffened  on  the  top  and 
bottom  sides 
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Figure  5-5.  Shear  stress  profiles  of  homogeneous  the  FGMs 
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Figure  5-6.  The  maximum  shear  stress  of  various  composites 


75  mm 
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Figure  5-7.  The  DSFGM  specimen.  A)  Without  matrix  cracks.  B)  With  matrix  cracks.  C) 
FE  model  to  demonstrate  the  SBS  test. 
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■  w/o  matrix  crack 

■  with  matrix  crack 


Figure  5-8.  Comparison  of  peak  loads  of  various  specimens  with  and  without  matrix 
cracks 


Figure  5-9.  The  GFC.  A)  Load-displacement  curve.  B)  Delamination  behavior  of  the 
specimen  without  matrix  crack  (d=6mm).  C)  Delamination  behavior  of  the 
specimen  with  matrix  crack  (d=6mm),  where  d  denotes  the  displacement  of 
the  indenter. 


Figure  5-10.  The  CFG.  A)  Load-displacement  curve.  B)  Delamination  behavior  of  the 
specimen  without  matrix  crack  (d=4mm).  C)  Delamination  behavior  of  the 
specimen  with  matrix  crack  (d=4mm). 
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Figure  5-1 1 .  The  LSQHybrid  composite.  A)  Load-dispiacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Deiamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 


Figure  5-12.  The  LSISFlybrid  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Deiamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 


Figure  5-13.  The  LS36Flybrid  composite.  A)  Load-displacement  curve.  B)  Deiamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Deiamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 
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Figure  5-14.  The  LSFGM  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Delamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 


Figure  5-15.  The  FS9Hybrid  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Delamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 


Figure  5-16.  The  FSISHybrid  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Delamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 
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Figure  5-17.  The  FS36Hybrid  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Delamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 


Figure  5-18.  The  FSFGM  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Delamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 


Figure  5-19.  The  DS9Hybrid  composite.  A)  Load-displacement  curve.  B)  Delamination 
behavior  of  the  specimen  without  matrix  crack  (d=6mm).  C)  Delamination 
behavior  of  the  specimen  with  matrix  crack  (d=6mm). 
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Figure  5-20.  Load-displacement  curves.  A)  DS18hybrid.  B)  DS36hybrid.  C)  DSFGM: 

The  locations  of  deiamination  in  the  symmetric  specimens  are  same  as  those 
of  the  DSOhybrid  shown  in  Figure  5-19. 
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Table  5-1 .  Various  hybrid  composite  laminates  and  FGM 


Name  Description 


Glass  fiber  composite 

Consists  of  unidirectional  glass  fiber  composites  (GFC). 

Carbon  fiber  composite 

Consists  of  unidirectional  carbon  fiber  composites  (CFC). 

LS9Hybrid  composite 

Consists  of  9  percent  CFC  and  91  percent  GFC.  CFC  was 
placed  on  the  loading  side. 

LSISHybrid  composite 

Consists  of  18  percent  CFC  and  82  percent  GFC.  CFC 
was  placed  on  the  loading  side. 

LS36Hybrid  composite 

Consists  of  36  percent  CFC  and  64  percent  GFC.  CFC 
was  placed  on  the  loading  side. 

FS9Hybrid  composite 

Consists  of  9  percent  CFC  and  91  percent  GFC.  CFC  was 
placed  on  the  free  side  (opposite  to  the  loading  side). 
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Table  5-1 .  Continued 


Name  Description 


FSISHybrid  composite 


FS36Hybrid  composite 


Consists  of  18  percent  CFC  and  82  percent  GFC.  CFC 
was  placed  on  the  loading  side  (opposite  to  the  loading 
side). 

Consists  of  36  percent  CFC  and  64  percent  GFC.  CFC 
was  placed  on  the  loading  side  (opposite  to  the  loading 
side). 


DS9Hybrid  composite 


Consists  of  9  percent  CFC  and  91  percent  GFC.  CFC  was 
placed  on  the  top  and  bottom  faces  (symmetric). 


DS18Hybrid  composite 


DS36Hybrid  composite 


Consists  of  18  percent  CFC  and  82  percent  GFC.  CFC 
was  placed  on  the  top  and  bottom  faces  (symmetric). 

Consists  of  36  percent  CFC  and  64  percent  GFC.  CFC 
was  placed  on  the  top  and  bottom  faces  (symmetric). 


LSFGM  composite 

Volume  fraction  of  CFC  is  varied  from  100  percent  at  the 
loading  side  to  0  percent  at  the  free  side  in  1 1  steps. 
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Table  5-1 .  Continued 


Name  Description 


FSFGM  composite  The  volume  fraction  of  CFG  at  layers  is  varied  from  100 

percent  at  the  free  side  to  0  percent  at  the  loading  side  in 
11  steps. 

DSFGM  composite  The  volume  fraction  of  CFG  is  varied  from  100  percent 

both  at  the  loading  and  free  side  to  0  percent  at  the  center 
in  11  steps  (symmetric). 


Table  5-2.  Elastic  properties  for  composite  materials  in  the  FEM 

E\=E2 

(GPa) 

E3 

(GPa) 

Gi3-G23=  Gi2 
(GPa) 

V13-V23 

V12 

Glass  fiber  composite 

26.9 

8.6 

3.1 

0.28 

0.16 

Carbon  fiber  composite 

70.0 

6.8 

3.3 

0.15 

0.13 
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CHAPTER  6 

CONCLUSIONS  AND  FUTURE  WORK 

Conclusions 

It  was  shown  that  the  damage  resistance  and  tolerance  of  laminated  composites 
can  be  enhanced  by  the  employment  of  translaminar  reinforcements  (TLR)  such  as 
stitching,  z-pinning  and  3D  weaving  and  also  by  hybrid  composites. 

A  non-dimensionai  anaiytical  model  focused  on  Mode  I  delamination  was 
developed  to  understand  the  role  of  the  TLR  on  delamination  behavior.  An  expiicit 
formula  for  the  apparent  interlaminar  fracture  toughness  was  derived  in  terms  of  the 
inherent  fracture  toughness  of  original  materials  and  the  bridging  force  due  to  z-pins. 
This  modei  is  capabie  of  estimating  the  apparent  fracture  toughness,  the  bridging  length 
and  allowable  bridging  force  thus  can  be  usefui  in  the  design  of  TLR  for  composite 
laminates. 

Along  with  understanding  advantages  of  TLR  in  increasing  the  damage  toierance 
of  iaminated  composites,  the  damage  behavior  of  iaminated  composites  subjected  to 
low  velocity  impact  ioading  was  studied.  Based  on  the  simiiarity  in  damage  deveiopment 
between  quasi-static  and  dynamic  loadings  observed  through  the  short  beam  shear 
(SBS)  tests,  the  FE  analyses  of  the  SBS  specimens  for  quasi-static  indentation  and  at 
severai  rates  of  iow-velocity  impact  ioadings  were  performed.  The  results  reveal  that 
inertia  effects  in  the  typical  velocity  range  of  the  striker  in  SHPB,  around  10m/s,  can  be 
negiigibie,  and  hence  the  quasi-static  anaiysis  is  usefui  and  valid  in  the  study  of  damage 
in  composite  specimens  under  iow  veiocity  impact  loading. 

The  delamination  behavior  of  3D  woven  composites  was  investigated  focusing  on 
the  effect  of  z-yarn.  The  3D  woven  composites  containing  both  single  and  double  z- 
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yarns  were  chosen  and  compared  with  the  2D  plain  woven  laminate.  The  double  z-yarn 
woven  composite  exhibited  enhanced  damage  tolerance  compared  to  the  single  z-yarn 
and  the  plain  woven  laminate.  The  relative  sliding  motion  between  two  layers  is 
constrained  by  z-yarns  thus  the  crack  propagation  of  delamination  is  suppressed.  This 
mechanism  increases  the  apparent  interlaminar  fracture  toughness  of  the  composites. 

The  interlaminar  shear  stress  profiles  in  the  various  hybrid  composites  were 
obtained  from  linear  FE  analyses  and  compared  within  the  framework  of  damage 
resistance.  Some  hybrid  configurations  resulted  in  reduced  maximum  shear  stress 
value  for  a  given  contact  force  thus  demonstrating  higher  damage  resistance.  Although 
matrix  cracks  may  change  the  location  of  delamination  and  decrease  damage 
resistance,  hybrid  composites  can  still  be  superior  in  structural  applications. 

Recommendations  for  Future  Work 

It  is  expected  that  more  studies  in  these  areas  would  provide  a  good  design  tool 
and  insights  into  application  of  composite  materials  in  structures  subjected  to  impact 
loading. 

First  of  all,  methodologies  for  accurate  prediction  of  intralaminar  failure  in  yarns 
and  matrix  regions  are  needed  and  their  implementation  in  the  FE  models  should  be 
studied.  Some  useful  experimental  data  provided  several  assumptions  for  using 
cohesive  elements  in  this  study.  However,  an  analysis  scheme  that  does  not  require 
a  priori  knowledge  of  crack  location  would  be  desirable. 

In  addition  to  matrix  cracks  and  delamination,  debonding  of  yarns  should  be 
considered  although  this  effect  is  less  than  delamination  in  orthogonal  woven 
composites.  The  study  including  debonding  can  be  used  in  modeling  of  different  types 
of  materials  such  as  angle-interlock  composites. 
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Finally,  hybrid  woven  composites  should  be  analyzed.  Each  layer  and  individual 
yarns  can  be  made  of  various  materials.  The  synergistic  effect  due  to  z-yarns  and 
hybridization  should  be  evaluated  for  the  best  performance  of  composite  materials 
subjected  to  impact  loading. 
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